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ABSTRACT 

We present high angular resolution Submillimeter Array (SMA) and Com- 
bined Array for Research in Millimeter-wave Astronomy (CARMA) observations 
of two GLIMPSE Extended Green Objects (EGOs)-massive young stellar object 
(MYSO) outflow candidates identified based on their extended 4.5 /xm emission in 
Spitzer images. The mm observations reveal bipolar molecular outflows, traced 
by high- velocity ^^C0(2-l) and HCO"'"(1-0) emission, coincident with the 4.5 fim 
lobes in both sources. SiO(2-l) emission confirms that the extended 4.5 /im emis- 
sion traces active outflows. A single dominant outflow is identified in each EGO, 
with tentative evidence for multiple flows in one source (Gil. 92— 0.61). The 
outflow driving sources are compact millimeter continuum cores, which exhibit 
hot-core spectral line emission and are associated with 6.7 GHz Class II CH3OH 
masers. Gil. 92— 0.61 is associated with at least three compact cores: the outflow 
driving source, and two cores that are largely devoid of line emission. In con- 
trast, G19.01— 0.03 appears as a single MYSO. The difference in multiplicity, the 
comparative weakness of its hot core emission, and the dominance of its extended 
envelope of molecular gas all suggest that G19.01— 0.03 may be in an earlier evo- 
lutionary stage than Gil. 92— 0.61. Modeling of the G19.01— 0.03 spectral energy 
distribution suggests that a central (proto)star (M ~10 M©) has formed in the 
compact mm core {Mgas ~12-16M0), and that accretion is ongoing at a rate of 
~10^^ Mq year^^. Our observations confirm that these EGOs are young MYSOs 
driving massive bipolar molecular outflows, and demonstrate that considerable 
chemical and evolutionary diversity are present within the EGO sample. 

Subject headings: ISM:jets and outflows — stars: formation — techniques: in- 
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Introduction 



Massive star formation remains a poorly understood phenomenon, largely due to the 
difficulty of identifying and studying massive young stellar objects (MYSOs)lj in the crucial 
early active accretion and outflow phase. During the earliest stages of their evolution, young 
MYSOs remain deeply embedded in their natal clouds. Most massive star-forming regions 
are also distant (> Ikpc) and crowded, with massive stars forming in close proximity to other 
MYSOs and to large numbers of lower-mass YSOs. Studying the early stages of massive star 
formation thus requires high angular resolution observations (to resolve individual objects 
in crowded regions) at long wavelengths unaffected by extinction. 

Large-scale Spitzer surveys of the Galactic Plane have yielded a promising new sample 
of young MYSOs with active outflows, which may be inferred to be actively accreting. Iden- 
tified based on their extended 4.5 //m er nission in Spitzer images , these sources are known 
as "Extended Green O bjects (EGOs)" JCvganowski et all boosl . boooh or "green fuzzies" 
(jChambers et al.ll2009l ) from the common coding of the 4.5 /im band as green in 3-color IRAC 
images. In active protostellar outflows, the Spitzer 4.5 fim broadband flux can be dominated 
by emission from shock-excited molecular lines (predominantly H2: Smith fc Rosen 20051 : 



De Buizer fc Vacca 



Smith et al.lbood : IPavis et al.lboOTi : lYbarra fc Ladalboool : lYbarra et al.lbolO 
20101). The resolution of Spitzer at 4.5 /im (~2") is sufficient to resolve the extended emission 
from outflows in massive star forming regions nearer than ~7 kpc. Over 300 EGOs have been 
cataloged in t he Galactic Legacy Infra red Mid-Plane Survey Extraordinaire (GLIMPSE-I) 
survey area by lCyganowski et al.l (120081 ). The mid-infrared (MIR) colors o f EGOs are consis- 
tent with those of young protostars still embedded in infalling envelopes (jCyganowski et al. 
20081 ) . A majority of EGOs are also associated with infrared dark clouds (IRDCs), identifled 
by recent studies as s i tes of the earliest stag es of massive star and cluster formation (e.g. 
Rathborne et allboO?! : Ichambers et al.lboOQh . 



Remarkably high detection rates for two diagnostic types of CH3OH masers in high- 
resolution Very Large Array (VLA) surveys provi de strong evidence that GLIMPSE EGOs 
are indeed massive YSOs with active outflows ( ICyganowski et al.l l2009l ). There are two 
Classes of CH3OH masers, both associated with star formation, but excited under differ- 
ent conditions by different mechanisms. Cla ss II 6.7 GHz CH 3OH masers are radiatively 

and references therein) and 



pumped by IR emission from warm dust (e.g 
are associated exclusively with ma ssive YSOs (e.g. 



Cragg et al.l 



2005 



Minier et al. 2003; Bourke et al. 2005 



Xu et al.ll2008l : iPandian et al.ll2008l ). Class I 44 GHz CH3OH masers are collisionally excited 



^We define MYSOs as young stellar objects (YSOs) that will become main sequence stars of M>8 Mq 
(O or early-B type ZAMS stars). 
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in molecular outflo ws, and particularly at interfaces between ou tflows and the surrounding 



ambient cloud (e.g. IPlambeck fc Mentenlll990l : Kurtz et al.l 12004 ). Of a sample of 28 EGOs 



>64% have 6.7 GHz Class II CH3OH masers (nearly double the detection rate of surveys 
using other MY SO selection criteria), an d of these 6.7 GHz maser sources, ~89% also have 
44 GHz masers (jCyganowski et al.l 120091 ). 



A complementary James Clerk Maxwell Telescope (JCMT; resolution ~20") molec- 
ular line survey towards EGOs with 6.7 GHz CH3OH maser detections found SiO(5-4) 
emissi on and HCO"'"(3-2) line proflles consistent with the presence of active molecular out- 



flows (jCyganowski et al.ll2009l ). SiO is particularly well-suited to tracing active outflows, 
as it persists in the gas pha se for only ~10^ years after being released by shocks (e.g. 
Pineau des Forets et al.l 119971 ). A single-dish (resoluti on ~80^0 3 mm sp ectral line survey 



of all EGOs visible from the northern hemisphere by IChen et al.l (I2OIOI ) found associated 
gas/dust clumps of mass 69-29000 M©, consistent with the identiflcation of EGOs as MYSOs. 
The nature of the driving sources of the 4.5 /im outflows is only loosely constrained by the 



^ „^ „ 

Cveanowski et al.l 


2009'). 


Cveanowski et al. 


(2009.) 



ering sources by the lack of VLA 44 GHz continuum detections 
A high detection rate (83%) for thermal CH3OH emission in the 
JCMT survey indicates the presence of warm dense gas, and possible hot core line emission. 

Further understanding of the nature of EGOs, and their implications for the mode(s) of 
high-mass star formation, requires identifying the driving source (s) and characterizing their 
physical properties, as well as those of the outflows associated with EGOs. Interferometric 
millimeter-wavelength line and continuum observations provide access to direct tracers of 
molecular outflows and dense, compact gas and dust cores, including a wealth of chemical 
diagnostics. In this paper, we present Submillimeter Array (SMA)g and Combined Array 
for Research in Mil limeter-wave Astr o nomy (CARMA)§ observations at 1 and 3 mm of 
two EGOs from the ICveanowski et al.l J2OO9I ) sample: Gl 1.92-0.61 and G19.01-0.03. The 
targets were chosen to have bipolar (and in some cases quadrupolar) 4.5 fim morphology, 
associated 24 fim emission, associated (sub) mm continuum emission in single-dish surveys. 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 

•^Support for CARMA construction was derived from the Gordon and Betty Moore Foundation, the 
Kenneth T. and Eileen L. Norris Foundation, the James S. McDonnell Foundation, the Associates of the 
California Institute of Technology, the University of Chicago, the states of California, Illinois, and Maryland, 
and the National Science Foundation. Ongoing CARMA development and operations are supported by the 
National Science Foundation under a cooperative agreement, and by the CARMA partner universities. 
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and 6.7 Class II and 44 GHz Class I CH3OH maser detections in the lCyganowski et al.l (|2009[ ) 
survey. The promise of extended 4.5 fim emission as a MYSO diagnostic lies largely in its 
ability to identify very young sources with ongoing accretion and outflow that are missed 
by other sample selection methods. These sourc es had not been ta r geted for study prior to 
their identification as EGOs and inclusion in the Cyganowski et al. ( 2009 1) sample, and very 
little is known about them beyond the results of that survey (see also §3.1.11 and §3.2.ip . 
In ^ we describe our observations, and in ^ we present our results. In §4] we discuss the 
physical properties of the compact cores and outflows associated with our target EGOs, and 
in §5] we summarize our conclusions. 



2. Observations 

2.1. Submillimeter Array (SMA) 

SMA observations of our target EGOs were obtained on 23 June 2008 with eight anten- 
nas in the compact-north configuration. The observational parameters, including calibrators, 
are summarized in Tabled! Two pointings were observed in a single track: Gil. 92— 0.61 
at a =18^13"^58M, 6 = -18°54'16'.'7, and G19.01-0.03 at a=18^25"^44^8, 6 = -12°22'45'.'8 
(J2000). The average 225 GHz opacity during the observations was ~0.26, with typical 
system temperatures at source transit of 220 K. In the compact-north configuration, the 
array is insensitive to smooth structures larger than ~20". The projected baseline lengths 
ranged from 7 to 88 kA. The double-sideband SIS receivers were tuned to a local oscillator 
frequency of 225.11 GHz, providing coverage of 219.1-221.1 GHz in the lower sideband (LSB) 
and 229.1-231.1 GHz in the upper sideband (USB). The spectral lines detected are reported 
in gSXland gSXl 

Initial calibration of the data was performed in MIRIAD. Each sideband was reduced 
independently, and the calibrated data were exported to AIPS. The AIPS task UVLSF was 
used to separate the line and continuum emission, using only line-free channels to estimate the 
continuum. The continuum data were then self-calibrated, and the solutions transferred to 
the line data. After self-calibration, the continuum data from the lower and upper sidebands 
were combined. Imaging was performed in CASA using Briggs weighting and a robust 
parameter of 0.5. The synthesized beam size is 3'.'2xl'.'8 (P.A.=59°) for G11.92 and 3'.'2xl'.'7 
(P.A.=63°) for G19.01. The la rms noise level in the continuum images is 3.5 mJy beam~^. 
The correlator was configured to provide a uniform spectral resolution of 0.8125 MHz. The 
line data were resampled to a velocity resolution of 1.1 km s~^, then Hanning smoothed. The 
typical noise level in a single channel of the Hanning-smoothed spectral line images is 100 
mJy beam~^. The ^^CO data were further smoothed to a resolution of 3.3 km s~^; the noise 
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in a single channel is 60 mJy beam ^. All measurements were made from images corrected 
for the primary beam response. 

Flux calibration was based on observations of Uranus and a model of its brightness 
distribution using the MIRIAD task smaflux. Comparison of the derived fluxes of the ob- 
served quasars (including 3C279, included as an alternate bandpass calibrator) with SMA 
flux monitoring suggests that the absolute flux calibration is good to <15%. The absolute 
position uncertainty is estimated to be 0'.'3. 

2.2. Combined Array for Research in Millimeter-wave Astronomy (CARMA) 

Our 3 mm CARMA observations were obtained on 29 July 2008 (Gil. 92-0.61) and 
30 July 2008 (G19.01— 0.03) in the D-configuration with 15 antennas (six 10.4 m and nine 
6.1 m antennas). The observational parameters, including calibrators, are summarized in 
Tabled! The projected baselines ranged from 1.5 to 31 kA for the 29 July observations and 
from 1.5 to 36.5 kA for the 30 July observations. The correlator was configured to cover SO 
(22-li) at 86.094 GHz (E„pper=19.3 K) in LSB, SiO (2-1) at 86.846 GHz (E„pper=6.3 K) in 
LSB and HCO+(1-0) at 89.189 GHz (E„pper=4.3 K) in USB with 31 MHz windows. Each 31 
MHz window consisted of 63 channels, providing a spectral resolution of 0.488 MHz (~1.7 
km s~^) and velocity coverage of ~100 km s~^. In addition, the correlator setup included 
two 500 MHz (pseudo) continuum bands, each comprised of 15 channels: one in LSB centered 
at ~85.7 GHz and one in USB centered at ~90.3 GHz. In the D-configuration, the array 
is insensitive to smooth structures larger than ~50". During the observations, the 230 GHz 
opacity ranged from ~0.47 to 0.5 on 29 July and from ~0. 46-0. 54 on 30 July. Typical 
(SSB) system temperatures at source transit were ~230-280 K on 29 July and ~190-240 K 
on 30 July. The phase center was a =18^13'^58M, 5 = -18°54'16'.'7 (J2000) on 29 July 
(Gll.92-0.61), and a =18^25"^44^8, 5 = -12°22'45'.'8 (J2000) on 30 July (G19.01-0.03). 

The data were calibrated in MIRIAD and imaged in CASA, using Briggs weighting 
and a robust parameter of 0.5. The synthesized beamsize is 6'.'8x5'.'l (P.A.=11°) for Gil. 92 
and 5'.'7x5'.'l (P.A.=-27°) for G19.01. Each band was reduced independently. The 3 mm 
continuum data for Gil. 92— 0.61 are not presented here because of aliasing from IRAS 
18110-1854, which is a much brighter continuum source than the target EGO. The IRAS 
source is ~1' northeast of the EGO, within the primary beam FWHP of the 6.1 m antennas 
(~130") but only at the 20% response level of the primary beam of the 10.4 m antennas 
(FWHP ~80"). For G19.01— 0.03, the continuum data from the upper and lower sidebands 
were combined to make a final continuum image with a la rms noise level of 0.5 mJy beam~^. 
All line data were Hanning smoothed to improve the signal to noise. The noise level in a single 
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channel of the Hanning-smoothed spectral line images is ~12 mJy beam""^ for Gil. 92— 0.61 
and ~9 mJy beam~^ for G19.01— 0.03. No continuum subtraction was performed, as the 
continuum contribution to each channel of the line data is negligible. For G19.01— 0.03, the 
3.4 mm continuum peak intensity (9.4 mJy beam~^) is less than the la rms level in the 
line datacube (10 mJy beam~^). For Gil. 92— 0.61, extrapolating the 1.3 mm continuum 
peak intensity to 3.4 mm (assuming a spectral index of three) likewise predicts a continuum 
contribution to the line data at the <la level. All measurements were made from images 
corrected for the convolved primary beam response of the heterogeneous CARMA array; in 
CASA this calculation is done in the visibility domain. 

Flux calibration was based on observations of Uranus and a model of its brightness dis- 
tribution using the MIRIAD task smafiux. Comparison of the derived fluxes of the observed 
quasars (including 3C273, observed as an alternate bandpass calibrator) with CARMA flux 
monitoring suggests that the absolute flux cahbration is good to <15%. 

In addition to the 3 mm observations described above, we obtained 1 mm observations 
of Gll.92-0.61 with CARMA in the C-conflguration on 25 Aprfl 2008 with eleven antennas 
(three 10.4 m antennas and eight 6.1 m antennas). The projected baselines ranged from 
9 to 184.5 kA, and the phase center was a =18^13'"58M S = -18°54'16'/7 (J2000). In 

the C-configuration, the array is insensitive to smooth structures larger than ~15". The 
(SSB) system temperature ranged from ~400-600 K during the observations. The correlator 
was configured to cover SiO(5-4) at 217.105 GHz and DCN(3-2) at 217.239 GHz in LSB 
and 80(56-45) at 219.949 GHz and CH30H(8o,8-7i,6) at 220.078 GHz in USB with 62 MHz 
windows and two 500 MHz (pseudo) continuum windows centered at ~216.45 GHz (LSB) and 
220.75 GHz (USB). Due to the high system temperatures and limited integration time, only 
the 500 MHz bands have sufficient signal-to-noise. The USB 500 MHz window encompassed 
the CH3CN(J=12-11) ladder, thus only the LSB (~216.5 GHz) (psuedo) continuum data are 
presented here. The absolute flux scale was set assuming a flux of 2.69 Jy for J1733-130 
based on quasar flux monitoring. The uncertainty in the absolute amplitude calibration is 
estimated to be ~20%. The data were calibrated in MlRlAD and imaged in CASA, using 
Briggs weighting and a robust parameter of 0.5. The resulting 1 mm continuum image has a 
synthesized beamsize of l'.'44x0'.'87 (P. A. =25°) and a la rms noise level of 4.3 mJy beam~^. 
As with the 3 mm data, all measurements were made from images corrected for the response 
of the heterogeneous primary beams. 
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Results 



3.1. Gll.92-0.61 



3.1.1. Previous Observations: Gil. 92— 0.61 



The EGO Gl 1.92—0.61 is ~1' SE of the more evolved massive star forming region 
IRAS 18110-1854. Single dish (sub)mm continuum maps t argeting the IRAS source show a 



millimeter /submillimeter clu mp coincident with the EGO ( iWalsh et al.ll2003l : iFaundez et al. 



2004J : [Thompson et al.ll2006l ). Strong (>240 Jy) H2O maser emission associated with the 
EGO was likewise serendipi tously detected in VLA observations targeting the IRAS source 
( iHofner &: Churchwelllll996l). The H2O maser source was s ubsequently included in the large 



single-dish CO survey of 



Shepherd fc Churchwelll (1l996l ). who detected broad ^^CO line 



wmgs. 



The MIR emission of Gil. 92— 0.61 is characterized by a bipolar 4.5 p,m morphol ogy 



with a NE and a SW lobe. The EGO is located in an IRDC fICvganowski et all I2OO8 



see 



also Fig. d^). The SW lobe is coincident with strong, blueshifted 44 GHz Class I CH3OH 
masers, while the NE lobe is sur rounded by an arc of systemic to slightly redshifted 44 GHz 
masers ( jCyganowski et al.ll2009l ). Elongated 24 fim emission is coincident with the NE 4.5 
lj,m lobe, as are tw o 6.7 G Hz Class H CH3OH masers (Fig. [T]). This EGO is unique among the 
Cyganowski et al.l ( 120091 ) sample in having multiple spatially and kinematically distinct loci 
of 6.7 GHz Class H CH3OH maser emission. The extended 24 /im morphology and multiple 
6.7 GHz maser spots suggest the possible presence of multiple MYSOs. SiO(5-4) and thermal 
CH30H(52,3-4i 3) emission were detected towards Gil. 92— 0.61 in s ingle-pointing JCMT ob - 
servations (resolution ~20") targeted at the NE lobe/24 /xm source (jCyganowski et al.ll2009l ). 
No 44 GHz continuum emission was detected towards the EG O to a 5cr sensitivity limit of 7 
mJy beam~^ ( resolution 0^.^99 x 0M4.ICvganowski et al.ll2009l ). We adopt the near kinematic 



distance from Cyganowski et al. ( 20091 ) for Gil. 92— 0.61: 3.8 kpc 



3.1.2. Continuum Emission: Gil. 92— 0.61 



Our 1.3 mm SMA and 1.4 mm CARMA data resolve three distinct compact contin- 
uum sources (Fig. [T^,b). These sources are designated MMl, MM2, and MM3 in order 
of descending peak intensity. Table [2] lists the observed properties of each source, includ- 
ing the integrated flux density and deconvolved source size determined from a single two- 
dimensional Gaussian fit. The (sub)mm clump coincident wit h the EGO is visible in the 
1.2 mm SEST/SIMBA map (resolution 24") of iFaundez et al.l ( 120041 ). No parameters are 
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tabulated, however, so we cannot compute the fraction of the single-dish flux recovered by 
the interferometers o 

All three mm continuum sources are coincident with the NE 4.5 /xm lobe. As shown in 
Figure [T]d, the MIPS 24 fim emission associated with the EGO is elongated along a N-S axis, 
and encompasses both MMl and MM3. The MIPS 24 /im image is saturated, introducing 
con siderable (^2-4") uncert ainty into the determination of the 24 fim centroid position (see 
also lCyganowski et al.ll2009l ). Figure [T]d shows, however, that the (saturated) 24/im peak lies 
~l-2" N of MMl, and roughly between MMl and MM3. The MIPS 24 /im counterpart thus 
likely consists of blended emission from these two s ources. MMl and MM3 a re also coincident 
with 6.7 GHz Class II CH3OH masers reported bv lCvganowski et all J2009h (Fig. [TJi,b). The 
CARMA 1.4 mm centroid position of MMl is offset ~0'.'6 (> 2100 AU) to the north of the 
intensity-weighted position of the southern 6.7 GHz CH3OH maser group (Gil. 918-0. 613), 
and the CARMA 1.4 mm centroid position for MM3 is offset ~0'.'4 (> 1500 AU) to the 
northeast of the intensity-weighted position o f the northern 6.7 GHz CH .sOH maser group 
(Gil. 919— 0.613). The H2O maser reported by lHofner fc Churchwelll (119961 ) is also coincident 
with MMl to within the astrometric uncertainties (Fig. Wp)- Notably, MM2 is offset to the 
northwest by >4" (> 14800 AU) from the 24 yum peak, and is not associated with a 6.7 GHz 
CH3OH maser. 



3.1.3. Compact Molecular Line Emission: Gil. 92— 0.61 

The continuum source MMl is associated with the richest molecular line emission in the 
Gil. 92— 0.61 complex. Emission in 27 lines of 11 species is detected towards Gil. 92— 0.61- 
MMl in our SMA observations. Table [3] lists the specific transitions, frequencies, and upper- 
state energies of lines detected at >3(T. Figure [2] shows the spectrum at the MMl continuum 
peak across the 4 GHz bandwidth observed with the SMA, with the transitions listed in 
Table [3] labeled. Table [3] also lists the peak line intensities, line centroid velocities, ^^fwhm^ 
and integrated line intensities obtained from single Gaussian fits to lines detected at >3o" at 
the MMl continuum peak. Some line profiles may be affected by outfiowing gas; transitions 
with non-Gaussian shapes are noted in Table [31 ^^CO and ^^CO were not fit, as the line 
profiles are complex and strongly self-absorbed. 



''Two major blind (sub)mm surveys of the Galactic Plane have recently been completed: the Bolo- 
cam Galactic Plane Survey (BGPS) at 1.1 mm and the ATLASG AL survey at 870 um. Unfortunately, 
Gll.92-0.61 falls outside the cover age of both the BGP S (lb| <0.5, iRosolowskv et al.lboiol ) and the 2007 
ATLASGAL campaign presented in ISchuUer et all (jiooi). 
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The spectrum of G11.92— 0.61-MMl is similar to those of hot cores observ ed with com- 
parable setups with the SMA, such as AFGL 5142 MM2 Jzhang et aliboOTh . The EGO 
spe ctrum is also similar to that of HH80-81 MMl, the driving source of the HH80-81 radio 
jet ( IQiu fc Zhang |2009| ). with the notable exception that 802(115,7-1248) emission (229.347 



GHz, E„pper=122 K) is not detected towards G11.91— 0.61-MMl. Emission from complex 
oxygen-rich organic molecules characteristic of strong hot cores (such as HCOOCH3) also is 
not detected towards Gil. 92— 0.61. 

Figures [3] and m present integrated intensity (moment zero) maps for selected transitions 
from Table O As shown in Figures [3] and HI emission from most species is compact and 
coincident with the continuum source MMl. Emission from high-excitation lines (E^pper^lOO 
K) is detected exclusively towards MMl. In contrast, the continuum source MM2 is devoid of 
line emission. The only species that exhibits compact emission coincident with the continuum 
source MM3 is C^^O (Fig. [3]). The CH3OH integrated intensity maps shown in Figure H] are 
discussed further in §3.1.51 

Most of the lines detected towards MMl are quite broad, with Avpwhm of 8-10 km s^^ 
The compact molecular line emission exhibits a velocity gradient, from SE (redshifted) to 
NW (blueshifted) (Fig. [5]). As shown in Figure [5], this gradient is consistent across species 
including SO, HNCO, CH3OH, and CH3CN. One possible explanation for this velocity gradi- 
ent is an unresolved disk, oriented roughly perpendicular to the outflow axis. Higher angular 
resolution data are required to investigate this possibility. Not all molecules detected towards 
the hot core show the same velocity gradient. One exception is OCS, which has redshifted 
emission to the NE and blueshifted emission to the SW. This is consistent with the kine- 
matics of the dominant outflow ( §3.1.41) . and suggests that the inner regions of the outflow 
may be contributing significantly to the observed compact OCS emission. 

Determining the vlsr^ of the mm continuum sources is complicated by the possibil- 
ity of confusion from outflowing gas or resolved-out emission from the extended envelope. 
For MMl, there is sufficient agreement among lines that exhibit compact emission and 
are detected with high signal-to-noise to estimate f2,s'/j(MMl)=35.2±0.4 km s^^. This is 
slightly blueward of the Vt.sr of 36 km s~^ estim ated from the lower angular resolution 
H^^CO^ observations of ICyganowski et al.l (120091 ). The systemic velocity of MMl is also 
blueshifted relative t o both the 6.7 GHz Cla ss II CH3OH masers coincident with MMl 
(v~37. 1 -37.6 km s~^: ICvea nowski et al.l 120091 ). and the peak H2O maser velocity (v= 



=40.7 



km s ^; iHofner fc Churchwelli 1 19961 ). There is also weaker H2O maser emission at the 6.7 



GHz CH3OH maser velocity. Table [3] lists a Gaussian fit to the C^^O emission towards the 
MM3 continuum peak. The emission is narrow (Avi?vKHM=3.5 km s~^), and has a line cen- 
troid velocity of 34.4 km s~^ Since no other compact line emission is detected associated 
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with MM3, however, it is difficult to be certain whether this velocity represents the MM3 gas 
vlsr- If so, then the thermal gas emission from MM3 is blueshifted by > 4 km s~^ relative to 
the coincident 6.7 GHz Class II CH3OH masers, which have velocities of ~38.6-39.5 km s^^. 
Since no emission centered on MM2 is detected, its vlsr cannot be determined. 



3.1.4. Extended Molecular Line Emission: Gil. 92— 0.61 



The observed low-excitation transitions of the abundant molecules ^^CO and HCO^ 
exhibit extended emission spanning a wide velocity range (> 80 km s~^) and most of the 
telescope field of view. Emission in SiO (2-1) is similarly spatially extended, but spans a 
narrower velocity range (~40 km s~^). While the kinematics of this extended emission are 
complex, the high velocity {\y—VLSR\ ^ 13 km s~^) ^^CO and HCO"^ emission are character- 
ized by a bipolar outflow centered on the continuum source MMl (Fig. [6]). To complement 
the integrated intensity maps of the high velocity gas shown in Figure 0, channel maps of 
the i2CO(2-l), HCO+(1-0), and SiO (2-1) emission are shown in Figures EE 

The red and blue lobes of the molecular outflow are asymmetric, both spatially and 
kinematically. The blueshifted lobe, SW of MMl, extends to more extreme velocities 
{\vmax,bine-VLSR\^5Q km s"^; \vmax,red-VLSR\^^Q km s"^) and further from the continuum 
source. The blueshifted lobe also exhibits stronger SiO (2-1) emission. The sense of the 
velocity gradient in the molecular gas agree s with that of 44 GHz C lass I CH3OH masers 
imaged with the VLA (resolution 0^^99xOM4. ICyganowski et al.ll2009l ). The concentration of 
blueshifted Class I masers coincides with the blueshifted molecular outflow lobe seen in ^^CO 
and HCO"'" and with the SW 4.5 lobe (Fig. Eti-jb). The Class I masers in the arc to the NE 
have near-systemic or slightly redshifted velocities, consistent with the location and more 
moderate velocity of the redshifted molecular outflow lobe. In particular, the SE section 
of the maser arc is coincident with moderately redshifted ^^CO, HCO"'", and SiO emission 
(Figs. [Tlini 39.9 and 46.5 km s~^ panels; this relatively low-velocity gas is not included in the 
integrated intensity maps shown in Figure E]) . 

The morphology and kinematics of the SiO (2-1) emission differ from those of the 
other outflow tracers, and copious SiO (2-1) emission is detected far from the mm con- 
tinuum sources (Figs. I6f9|) . The excitation of low- J rotational lines of SiO, such as the 
2-1 transiti on, depends priniarily on the density, n^^ ( as opposed to the kinetic tempera- 



ture, Tkin] iNisini et al.l 120071 : Ijimenez-Serra et al.ll2010l ). and extended (parsec-scale), qui- 



escent (Av~0.8 km s~^) S iO(2-l) emission has recently been observed towards an IRDC 
( iJimenez-Serra et al.ll2010l ). The comparatively broad linewidths of the SiO(2-l) emission 
towards Gil. 92— 0.61 indicate that the entirety of the observed SiO emission is attributable 
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to outflow-driven shocks, with bright SiO(2-l) knots hkely tracing the impact of these shocks 
on dense regions in the surrounding cloud. 

While the NE(red)-SW(blue) gradient dominates the ^^CO and HCO"*" velocity fields, 
there are other features that suggest multiple outflows may be present. In particular, 
blueshifted ^^CO, HCO+, and SiO emission are detected NE of MMl at velocities > 10 
km s~^, and redshifted emission SW of MMl at velocities < 65 km s"^ (Figs. |6]|9]). Near the 
Vlsr of ~35 km s~^, however, low- velocity outflow emission is confused with emission from 
ambient gas. 

Emission from 80(65-54) also extends NE and SW of MMl (Fig. [3]). The morphology 
and kinematics are consistent with this SO emission arising in the dominant outflow; in 
contrast to SiO, the SO emission is stronger towards the redshifted (NE) outflow lobe. 
The properties of the lower-excitation SO (22-I1) emission (E„pper=19.3 K) observed with 
CARMA (not shown) are similar to those seen in SO (65-54) at higher spatial and spectral 
resolution with the SMA. Faint SO (65-54) emission is detected coincident with the MM3 
continuum source, at a velocity consistent with that of the C^^^O ( §3.1.3p . However, since 
the bipolar outflow(s) overlap the MM3 position, it is unclear whether this SO emission is 
associated with MM3. 



3.1.5. Millimeter CH^OH masers: Gil. 92- 0.61 



As shown in FigureH the morphology of the 229.759 GHz CH30H(8_i,8-7o,7)E (E^ 



upper " 



=89 

K) line emission is strikingly different from that of any other observed CH3OH transition. 
There is very strong CH30H(8_i^8-7o,7) emission to the southwest of MMl, coincident with 
the brightest 44 GHz CH30H(7o-6i)A+ Class I masers detected bv lCvganowski et all J2009I ) 
(Fig. H]). Strong 229.759 GHz CH3OH emission is also observed NE of MMl, also coincident 
with 44 GHz Class I masers. The 44 GHz (7n-6i ) and 229.759 GHz (8_i,8-7o,7) lines are both 
Class I CH3OH maser transitions. Isivsh et al.l J2002h first reported 229.759 GHz CH3OH 
maser emission towards DR21 (OH) and DR21 West based on observations with the IRAM 
30 m telescope. Probable maser emission in this transitio n has been detected with the SMA 
in HH 80-81, coincident with a 4 4 GHz CH3OH maser Joiu fc Zhanelboogh . and in IRAS 
05345+3157 jFontani et al.ll2009h . 



The CH30H(8_i,8-7o,7) emission observed NE and SW of G11.92-0.61-MM1 is spec- 
trally narrow (Fig. Hb-c). In both cases, the velocity of the mm emission agrees well with the 
velocities of the coincident 44 GHz masers. The extended appearance of the northeastern 
229.759 GHz emission is consistent with emission from multiple masers being blended at the 
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lower spatial resolution of the SMA obser v ation s (^ 3^^ compared t o ^0^/5 for the 44 GHz 
VLA data). Like the studies of ISlysh et al.l ( 120021 ) and lQiu fc Zhana ( 120091 ). the beamsize of 
our observations is too large to definitively establish the maser nature of the emission based 
on its brightness temperature. The peak intensity of the SW CH30H(8_i^8-7o,7) emission 
is 3.5 Jy bea m~^, correspondin g to Tb=14.3 K. For the NE emission, Ipeafe=l-6 Jy beam~^ 
(Tb=G-G K). ISlysh et al.l ( 12002| ) use the 229.759/230.027 line ratio as a discrimant between 
thermal and masing 229.759 GHz emission, with ratios > 3 indicative of non-thermal 229.759 
GHz CH3OH emission. The 230.027 GHz line is a Class II transition, so is not expected to 



be in verted under conditions that excite 229.759 GHz Class I maser emission (jSlysh et al. 



2OO2I). The 229.759/230.027 ratios for the SW and NE spots in Gil. 92-0.61 are 100 and 7, 



respectively. For comparison, the ratio at the MMl 1.3 mm continuum peak is 2, consistent 
with thermal emission from the hot core. The SW and NE 229.759 GHz emission features 
coincide spatially and spectrally with 44 GHz Class I CH3OH masers; this agreement, along 
with the millimeter line ratios and narrow linewidths, strongly supports the interpretation 
of these emission features as mm CH3OH maser emission. 



3.2. G19.01-0.03 



3.2.1. Previous Observations: G19.01-0.03 



This source was entirely unknown prior to being cataloged as an EGO by lCyganowski et al. 



( 120081 ). The EGO G19.01— 0.03 has a striking MIR appearance, with bipol ar 4.5 Atm emission 



centered on a point source detected in GLIMPSE and MIPSGAL images. ICyganowski et al. 



( I2OO9I ) detected kinematically complex Class II 6.7 GHz CH3OH maser emission coincident 
with the "central" source. Copious 44 GHz Class I CH3OH maser emission is associated with 
the 4.5 /xm lobes, with blueshifted masers con centrated to the north of the MIR point source 
and systemic /redshifted masers to th e south (ICyganowski et al.l 120091 ). The EGO is located 
in an IRDC (ICyganowski et al.ll2008l ). and a (sub)mm continuum source coincident with the 
EGO is detected in blind single-dish surveys of the Galactic Plane at 1.1 mm and 870 /xm 
(BGPS and ATLASGAL, iRosolowskv et al.lboid ISchuUer et aLll2009h . SiO(5-4) and ther- 
mal CH30H(52^ 3-4i 3) emission were detected towards G19.01— 0.03 in single-pointing JCMT 
observations (resolution ~20") targeted at the MIR point source ( ICyganowski et al.ll2009l ). 
No 44 GHz continuum emission was detected towards the EG O to a 5cr sensitivity limit of 5 
mJy beam~^ ( resolution 0^^59 x ^ ^51.ICvganowski et al.ll2009l ). We adopt the near kinematic 



distance from ICveanowski et al.l (l2009|) for G19.01-0.03: 4.2 kpc. 
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3.2.2. Continuum Emission: G19.01-0.03 



At the resolution of our SMA 1.3 mm and CARMA 3.4 mm observations, the miUimeter 
continuum emission associated with the EGO G19.01— 0.03 appears to arise from a single 
source, called MMl. The SMA 1.3 mm continuum image is shown in FigureHt. The position 
of the peak 3.4 mm continuum emission (not shown) coincides with the 1.3 mm continuum 
peak. In our SMA 1.3 mm continuum image, we recover 7.6l}'4% of the single-dish flux 
density of 3.6^n' fi Jy measured from the 1.1 mm Bolocam Galactic Plane Survey (BGPS, 



resolution ~30^^ iRosolowskv et al.lboioh FI 



The fitted position of MMl from the SMA 1.3 mm continuum image is coincident with 
the i ntensity-weighted position of the Class II 6.7 GHz CH3OH maser ( ICyganowski et al. 
2009h and with the GLIMPSE point source SSTGLMC G019.0087-00.0293 within the ab- 
solute positional uncertainty of the mm data (0'.'3). The MIPS 24 /im peak is offset to the 
NW by ~1'.'3 (> 5500 AU). This is consistent, however, with MMl being coincident with 
the MIPS 24 /xm source w'ithin the absolute positional uncertainty of the MIPSGAL survey 
(median 0'.'85, up to ~3", ICarev et al.ll2009h . 



3.2.3. Compact Molecular Line Emission: C19.01-0.03 

Emission in 15 lines of 9 species is detected towards G19.01— 0.03-MMl in our SMA 
observations. Table H] lists the specific transitions, frequencies, and upper-state energies of 
lines detected at >3cr. Figure dO] shows the spectrum at the MMl continuum peak across 
the 4 GHz bandwidth observed with the SMA, with the transitions listed in Table H] labeled. 
As shown by a comparison of Figures |2] and [101 many of the stronger lines detected towards 
G11.92-0.61-MM1 are also detected towards G19.01-0.03-MM1. These lines are much 
weaker towards G19.01— 0.03, however, and emission from higher-energy transitions {Eupper> 
200 K) is notably lacking. The highest-energy line detected towards G19.01— 0.03-MMl is 
the k=4 component of the CH3CN(12-11) ladder (E„pper=183 K), and this line is very weak 
(~3.5cr). Other than the k=3 and k=4 CH3CN lines, no line emission with Eupper->100 K is 
detected towards G19. 01— 0.03-MMl. As shown in Figures ITT] and WI\ emission from most 
species is compact and coincident with the continuum source. 

Table m also lists the peak line intensities, line center velocities, ^'VpwHM, and integrated 
line intensities obtained from single Gaussian fits to lines detected at >3(j at the G19.01— 0.03 



^BQps fl^x dens ity from 
Dunham et"al] (|20ld ). 



Rosolowskv et al. 



(|2010( ) with correction factor apphed as discussed in 
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MMl continuum peak. Some line profiles may be affected by outflowing gas; lines not well 
fit by a Gaussian are noted in Table H As for Gll.92-0.61, the ^^CO and ^^CO lines 
were not fit, as the line profiles are non-Gaussian and, in the case of ^^CO, strongly self- 
absorbed. As Table S] demonstrates, there is good agreement among the central velocities 
determined from the different species and transitions. Considering all lines detected at >5cr, 
'yLSR(MMl)=59.9=bl.l km s~^. This is in good agreement with the central velocities of 
Hi3CO+(3-2) (vcenter=59.9±0.1 km s-i) a nd CH.,OH('59 .-^-4i Q ('v. ,„,,.=59.7±0.2 km s"^) 



observed with the JCMT (resolution ~20", ICyganowski et al.l 120091 ). 



Compared to G11.92-0.61-MM1, the lines detected towards G19.01-0.03-MM1 are 
relatively narrow. Most of the transitions detected with >5a have Avfwhm < 4 km s~^. The 
6.7 GHz Class II CII3OH masers associated with G19.01— 0. 03-MMl span a compar atively 



wide velocity range of ~7.5 km s~^, from 53.7-61.1 km s~^ (ICyganowski et al.ll2009l ). The 
velocity range of the 6.7 GHz CH3OH maser emission extends much further to the blue of 
the vlsr than it does to the red. 

Figures [TT] and [T2] present integrated intensity (moment zero) maps for selected tran- 
sitions in Table ID The only species in Figure [TT] that exhibits significant emission not 
coincident with the continuum source is C^^O (2-1). (The properties of the SO (22-I1) 
(Eupper=19.3 K) emission observed with CARMA (not shown) are similar to those of the 
SO (65-84) emission observed at higher spatial and spectral resolution with the SMA.) The 
C^^O emission coincident with the continuum source has near-systemic velocities, while the 
two knots of C^^O emission to the south of MMl are redshifted and likely associated with 
knots in the outflow. 



3.2.4- Extended Molecular Line Emission: G19.01—0.03 

Extended molecular line emission, spanning most of the telescope field of view, is exhib- 
ited by i2CO(2-l) and HCO+(1-0). Extended SiO(2-l) emission is also observed. Figure [H] 
presents integrated intensity images of high- velocity gas, while Figures [TU [T51 and HH] show 
channel maps of the ^^C0(2-l), HCO"'"(1-0), and SiO(2-l) emission, respectively. 

As shown in Figures [T3lfT5| ^^CO and HCO+ trace a bipolar molecular outflow centered 
on the continuum source MMl. The outflow axis is roughly N-S, with the blueshifted lobe to 
the north and the redshifted lobe to t he south. This is consist ent with the velocity gradient 



of the 44 GHz Class I CH3OH masers f ICvganowski et al.ll2009r ). Blueshifted 44 GHz CH3OH 



masers are concentrated towards the northern 4.5 /xm lobe, while systemic and redshifted 44 
GHz CH3OH masers are concentrated to the south of the central source. The 44 GHz Class 
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I CH3OH masers trace the edges of the high- velocity ^^CO and HCO"'" lobes remarkably well 
(Fig. [T3b.b). In particular, an arc of 44 GHz masers appears to trace the edges and terminus 
of the blueshifted ^^CO jet. 

The full velocity range of the ^^CO outflow is ~135 km s~^. The kinematics of the out- 
flow are notably asymmetric. The highest-velocity blueshifted gas has \vmax,biue—VLSR\^^OG 
km s^^, while for the redshifted lobe, \ymax,red—'VLSR\ is only ~29 km s~^. The veloc- 
ity distribution of the 44 GHz CH3OH masers is also asymmetric with re spect to the Vr,.<iR'- 

''Vmm:,blne,maser-VLSRhQ-5 km S"^ whllc \Vmax,red,maser-VLSRh'2-'2 km S"^ (jCygaUOWski et al. 



20091 ). The blueshifted lobe of the molecular outflow traced by ^^CO and HCO"*" is coincident 
with the northern lobe of extended 4.5 /xm emission. In contrast, the most highly redshifted 
molecular gas is found south of the brightest 4.5 fim emission in the southern lobe. The high 
velocity outflow gas is clumpy. Both the red and blue lobes are characterized by strings of 
bright knots. 

The SiO(2-l) emission differs in kinematics and morphology from the ^^CO and HCO^ 
emission (Fig. [T3]). Very little redshifted SiO emission is detected. Blueshifted SiO emis- 
sion is concentrated north of the continuum source, consistent with the orientation of the 
^^CO/HCO^ outflow. The morphology of the blueshifted SiO emission, however, is linearly 
extended along an E-W axis. Near the systemic velocity (~60 km s~^), the SiO emission 
extends to the south towards the continuum source (Fig. [T6l) . Like the high- velocity ^^CO 
and HCO"^ emission, the SiO emission is characterized by clumps and knots. The strongest 
blueshifted SiO emission arises from a clump offset to the east of the extended 4.5 /xm 
emission. 

Near the ul5_r(~60 km s~^), the ^^CO image cube shows artifacts from large-scale emis- 
sion resolved out by the interferometer (Fig. [T^ . This suggests that near the systemic 
velocity, the ^^C0(2-l) emission is dominated by emission from a large-scale extended en- 
velope. This interpretation is consistent with the large spatial extent of the HCO"'"(1-0) 
emission near the vlsr (Fig. [T5|l . 



3.2.5. Millimeter CH^OH masers: G 19. 01- 0.03 



As seen in Gll.92-0.61, the 229.759 GHz CH30H(8_i,8-7o,7)E (E„ppe^=89 K) emission 
towards G19.01— 0.03 has a very different morphology than any other observed CH3OH 
transition. Figure [T2] presents integrated intensity maps of the 230.027 GHz (E„pper=40 K), 
229.759 GHz, and 220.078 GHz (E„ppe^=97 K) CH3OH emission t oward s G19.01-0.03, with 
the positions of 44 GHz Class I masers from ICyganowski et al.l (120091 ) marked. There are 
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three compact loci of 229.759 GHz CH3OH emission to the north of MMl, two of which 
are coincident with numerous 44 GHz masers. Figure [12] also shows the profiles of the 
230.027 GHz, 229.759 GHz, and 220.078 GHz CH3OH emission towards the peak of each 
of these spots. The 229.759 GHz CH3OH emission towards these loci is spectrally narrow. 
Towards the two northern spots, the velocity of the 229.759 GHz CH3OH emission agrees 
well with the velocity range of the coincident 44 GHz CH3OH masers. In the integrated 
intensity map shown in Figure [12], the southernmost of the three 229.759 GHz CH3OH 
emission spots appears to lie between three clusters of 44 GHz CH3OH masers. The SMA 
data has much lower spatial and spectral resolution than the 44 GHz VLA data (SMA: 3", 
1.6 km s~^; VLA: 0'.'75, 0.24 km s~^). Examination of the data cubes shows that the 229.759 
GHz CH3OH emission near this location is consistent with being a blend of masers seen 
at 44 GHz, given the lower spatial and spectral resolution of the mm data. As shown in 
Figure [T^ the 229.759/230.027 line ratios to wards the three sp ots are all >3, consistent with 
non-thermal 229.759 GHz CH3OH emission ( ISlysh et al.ll2002l . see also §3.1.5p . Towards the 
MMl continuum peak, this line ratio is 2, consistent with thermal emission. 



Discussion 



4.1. Continuum Sources 



4- 1.1. Spectral Energy Distributions (SEDs) 

Unfortunately, the three members of the Gil. 92— 0.61 (proto) cluster are not resolved in 
existing data at wavelengths shorter than 1.3 mm. To better constrain the SED of the proto 



cluster as a whole, we measured the 70 /xm flux density from MIPSGAL images (jCarey et al. 



-50% due to 



(Carev et al. 


2009 


)• 


Walsh et al. 


(2003) 



port integrated flux densities of 140 Jy (450 /xm) and 12 Jy (850 /xm) for the (sub)mm 
clump (G11.92— . 64B in their nomenclature). Using these fluxes and the 24 /xm flux from 



f l2007h \i 



Cvganowski et al.l ( l2009l ). we fit the 24-850 /im data using the model fitter of iRobitaille et al. 



We do not include IRAC data in the SED because emission mechanisms not in- 
cluded in the models (e.g. line emission from shocked molecular gas and PAHs) may con- 
tribute significantly to the IRAC bands in this source. The fits, shown in Figure [T7k . are 
consistent with a bolometric luminosity o f ~10^ for the cluster as a whole (Fig. [T7b). 



The models in the iRobitaille et al.l (120071 ) grid assume a single central object in determing 



http://caravan.astro.wisc.edu/protostars/ 
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source properties (such as stellar mass) from the SED. Since Gil. 92— 0.61 is a (proto) cluster 
and unresolved at A < 1.3 mm, the model fitter cannot be used to constrain the properties 
of the cluster members. 

The EGO G19.01— 0.03 is unusual in that the "central" source is clearly resolved from 
the exten ded emission in IRAQ im ages (unlike most EGOs), and is a GLIMPSE point source 
(see also ICyganowski et al.ll2008l ). SED modeling can thus be used to infer the properties 
of the source driving the 4.5 /im outflow. As for Gil. 92— 0.61, we measured the 70 fim 
flux density of G19.01-0.03 from the MIPSGAL image. We find a flux density of 223 
Jy, with an estimated uncertainty of 50%. Using this measurement, GLIMPSE catalog 
photometry for the point source SSTGLMC G019.0087-00.0293, the 24 /xm flux density from 



Cvganowski et al.l (120091) . and the 1.3 mm flux density from Tabled we fit the SED using the 
Robitaille et all (12003) model fitter. As shown in Figures [TMTgt the SED is well-fit (Fig.UHJi) 



by models with bolometric luminosity of ~10^ Lq (Fig. [T8b). stellar mass ~10 M© (Fig. [T9|) . 
envelope mass ~10^ M0 (Fig. [T8b). and envelope accretion rate ~10~^ M0 yr~^ (Fig. [T9|) . 
Stellar age is also a pararneter o f the models, but is not well-constrained (Fig. [TSld). In the 
scheme of iRobitaille et al.l (120061 ) . evolutionary stage is defined by the ratio of the envelope 
accretion rate to the stellar mass. The youngest sources. Stage O/I, are defined as having 
Menv/M^ >10-^ yr-1. For G19.01-0.03, Menv/M^^IQ-"^ yr'^ (Fig. [19]), indicating that the 
source is likely to be very young. 



4- 1-2. Temperature Estimates from Line Emission 



The J =12-ll CH.sCN ladder is well-suited for measuring the ga s temperature in ho t 
cores (e.g. IPankonin etaPboOll : lArava et alJbood : Izhang et al]l2007l : lOiu fc Zhangl l2009[ ) . 
Figures [20] and [21] show the best-fit single-component model of the CH3CN emission overlaid 
on the observed spectrum at, respectively, the G11.92-0.61-MM1 and G19.01-0.03-MM1 
continuum peak. For each CH3CN emission component, the mode|l| assumes local ther- 
modynamic equilibrium (LTE) and the same excitation conditions for all K components, 
and accounts for optical depth effects and emission from the isotope CHg'^CN. The veloc- 
ity (frequency) separations of the K components are fixed to the laboratory values. The 
temperature, size (diameter), and CH3CN column density of the emitting region are free 
parameters, and the model that best fits the observed spectrum is found by minimizing the 
mean squared error. The parameters of the best-fit models are summarized in Table [3 

A single-component model provides an adequate fit to the CH3CN spectrum of G19.01— 0.03- 



^Developed using the XCLASS package, 



http://www.astro.uni-koeln.de/projects/schilke/XCLASS 
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MMl (Fig. [2T| T~114 K, size~2500 AU). In contrast, a single- component model is a no- 
tably poor fit to the CH3CN spectrum of G11.92-0.61-MM1 (Fig. [20]). In particular, the 
model severely underpredicts the emission from the k=7 and the (blended) k=0/l lines, 
while overpredicting the emission from most of the intermediate k components (k=3,4,6). 
To investigate this discrepancy, we allowed for two CHsCN-emitting regions, with different 
temperatures, sizes, and column densities. As shown in Figure [201 ^ two-component model 
that includes a compact (~2300 AU), warm (~166 K) component and an extended (~11400 
AU), cool (~77 K) component provides a much better fit to the data (see also Table \5^. 
This combination of parameters is likely not unique, and certainly we expect that the real 
emission exhibits a gradient in temperature rather than a step function. Even so, this result 
convincingly demonstrates that both cool and warm temperatures are present. Interestingly, 
the physical scale of the warm component (~2300 AU) agrees remarkably well with that of 
CHgCN-emitting region in G19.01-0.03-MM1 (~2500 AU). 

Five transitions of CH3OH are detected towards the G11.92— 0.61-MMl continuum peak, 
with Etjpper=40-579 K. This is sufficient to obtain an ind ependent estimate of the g as tem- 



perature by applying the rotation diagram method (e.g. [Goldsmith fc Langerl 119991 ) to the 
observed CH3OH emission, using the relations: 



Nu 3k 



Tdv, (1) 



and 

log{N^)/{gu) = log{Ntot/QiTrot)) - 0A3A3EjkTrou (2) 

where is the column density in the upper state, k is Boltzmann's constant, u is the line 
rest frequency, g/ is the nuclear spin degeneracy, gx is the K degeneracy, /i^S is the product 
of the square of the molecular dipole moment and the line strength, J Tdv is the observed 
integrated intensity of the line, Ntot is the total column density, Q(T rot) is the partition 
function evaluated at the rotation temperature Trot, and E„ is the upper state energy of the 
transition. For each transition, the integrated line intensity was determined from a single 
Gaussian fit to the line emission at the 1.3 mm continuum peak (Table [3]) • Both A and E 
transitions are included in the rotation diagram analysis; we do not detect enough transitions 
to treat the two types separately. The rotation temperature derived from a weighted least- 
squares fit to the data is 230±39 K (Fig. [22|) . somewhat higher than the temperature derived 
from the CH3CN fitting. 



As discussed in detail in [Goldsmith fc Langen ( 1l999[ ). however, optical depth effects can 



inflate the temperature de rived from a rotation diagram analysis. We follow the procedure of 
Brogan et al. J2007I . l2009h in iteratively solving for the values of r and Trot that best fit the 



data {Cr[i] = r[z]/(l — e"^'*') where i refers to the ith spectral line). As shown in Figure 
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this improves the fit considerably. The optical depth in the line with the highest opacity 
(CH30H(8_i,8-7o,7) at 229.759 GHz) is 3.67. With the optical depth correction, the derived 
temperature is 166±19 K. This is in remarkably good agreement with the temperature of the 
warm component from the CH3CN analysis (166 K). Only three CH3OH transitions (Eupper 
40-97 K) are detected towards G19.01— 0.03-MMl, too few for accurate rotation diagram 
analysis, but consistent with the cooler temperature derived for this source from CH3CN. 



4.1.3. Mass Estimates from the Dust Emission 

At millimeter wavelengths, thermal emission from dust and free-free emission from ion- 
ized gas can both contribute to the observed continuum flux. For our target EGOs, the 
available limits on any free-free contribution are not very stringent. Neither Gil. 92— 0.61 
nor G19.01— 0.03 had detectable 44 GHz continuum emission in the VLA observations of 



Cyganowski et al.l ( l2009l ). The 5a limits are 7 mJy beam~^ (synthesized beam 0'.'99x0'.'44) 
and 5 mJy beam~^ (synthesized beam 0'.'69x0'.'51) respectively. Extrapolating the 5a 44 
GHz upper limits assuming optically thin free-free emission [S^ oc z/", a=-0.1) gives upper 
limits at 1.3 mm of 6.0 mJy for Gil. 92-0.61 and 4.3 mJy for G19.01-0.03. For the adopted 
dust temperatures for G11.92-0.61-MM1 and G19.01-0.03-MM1, a free-free contribution 
at this level would have a minimal impact on the mass estimates (^0.4 Mq). If we instead 
extrapolate the 5a 44 GHz upper li mits assumin g a spectral index a=l (appropriate for a 



hypercompact (HC) HII region, e.g. iKurtzl 120051 ). the effect on the mass estimates is more 



substantial, up to ~2.5 Mq. For the weakest mm continuum source, G11.92— 0.61-MM3, 
free- free emission from a HC HII region (a=l) could in principle account for the entirety of 
the 1.4 mm flux density observed with CARMA and the majority (~73%) of the 1.3 mm 
flux density observed with the SMA. Deep, high-resolution continuum data at a range of 
cm wavelengths are required to constrain the presence and properties of any ionized gas 
associated with our target EGOs. In the absence of available evidence to the contrary, 
we assume the entirety of the millimeter-wavelength continuum emission is attributable to 
thermal emission from dust. 

Table [6] presents estimates from the thermal dust emission for the gas mass Mgas, column 
density of molecular hydrogen Nh2, and volume density of molecular hydrogen nH2, for 
G11.92-0.61-MM1, MM2, and MM3 and G19.01-0.03-MM1. The gas masses are calculated 
from: 

_ 4.79 X 10-''RSAJy)D\kpc)Cr,^^, 

where R is the gas-to-dust mass ratio (assumed to be 100), S,^ is the integrated flux density 
from Table [21 D is the distance to the source, C^-^^^^ is the correction factor for the dust 
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opacity C 



'^dust 



opacity coefficient in units of cm^ g ^. For gas densities of 10^-10^ cm ^, ^^^ ^>.mrn.^l for dust 



: dust ) is the Planck function, and is the dust mass 

grains with thick or thin ice mantles (lOssenkopf fc Henning||l994j ). Scaling from fi;i 3mm=l 
assuming /3=1.5, we adopt fi;3,4mm=0.24. We estimate a range of dust temperatures for each 
source based on its observed spectral hne properties (discussed in detail below). The dust 
opacity, Tdust = — ln(l — 7^^), is derived using the beam-averaged brightness temperature 

-^dust 

(Tfe) and assumed dust temperature (Tdust) for each source and listed in Table [61 The 
calculated values of T^ust are generally small (<0.1), indicating that the dust emission is 
optically thin. The column densities and volume densities presented in Table [6] are also 
beam-averaged values. 

As noted above, estimating gas masses using equation [3] requires an estimate of the 
dust temperature. For G11.92-0.61-MM1 and G19.01-0.03-MM1 we use the values of 
Tgas derived from the CH3CN and CH3OH emission ( §4.1.2p . At the high gas densities 
implied by our obser vations (>10^ cm~^), the gas and dust tern peratures are expected to 
be well-coupled (e.g. ICeccareUi et~aDll996l : iKaufman et"aDll998h . For G11.92-0.61-MM1, 
the situation is complicated by the presence of two temperature components, implied by 
the CH3CN fits ( §4.1.2p . Both the compact warm (size ~0'.'6) and more extended cool (size 
~3'.'0) emission regions are similar in scale to the 3'.'2 x I'.'S SMA beam. A step-function 
temperature structure is physically unrealistic, but the sensitivity and spatial resolution of 
the present observations are insufficient to resolve the temperature gradient in MMl. In the 
future, the sensitivity and high spatial resolution attainable with ALMA will allow detailed 
investigation of the temperature structure. Since the observed millimeter continuum is likely 
a mix of emission from the warm and cool components, we adopt a broad temperature range 
(70-190 K) for the estimates in Table [61 

Constraining the temperatures of G11.92-0.61-MM2 and G11.92-0.61-MM3 is more 
difficult because of the paucity of associated line emission. MM2 lacks clear MIR counterparts 
at 3.6-24 fim, is completely devoid of mm- wavelength line emission, and has no known 
maser emission. In contrast, MM3 emits at 24 /im and is associated with 6.7 GHz Class II 
CH3OH masers and possibly with a C^'^O clump. MM3 is also associated with the brightest 
8 yum emission in the region (Fig. [I^,b). Taken together, the evidence strongly suggests 
that MM3 is warmer than MM2. For MM2, we adopt a temperature range Trf„st=20-40 K 
based on the absence of associated molecular line emission. The 6.7 GHz CH3OH masers 
associated with MM3 are quite weak (peak T5~16500 K, 1'.'94 x ?96 synthesized beam), a s 
are the CH3OH masers associated with MMl (peak T5~7400 K, ICyganowski et al.ll2009l ). 
Class II CH3OH masers a re radiatively pump ed by infrared photons emitted by warm dust 
(e.g. 



Cragg et all 120051 ) 



Cragg et al.l ( 119921 ) found that a blackbody with T<50 K was 
sufficient to excite moderate 6.7 GHz CH3OH maser emission (T;, <6xl0^ K). More detailed 
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investigations of Class II CH3OH maser excitation have focused pri marily on the para meter 



space that gives rise to bright (T;, >10^ K) maser emission (e.g. ICragg et al.l l2005l who 



invoke dust temperatures >100 K). No high-excitation molecular lines (Eupper>100 K) are 
observed towards MM3. In sum, the multiwavelength data suggest two possibilities: MM3 
may be of intermediate temperature, or may be hotter (and more evolved) and simply have 
very little molecular material left around it. Additional data are required to constrain the 
nature and evolutionary state of MM3 ( §4.3.2p : we adopt a range of Trf„st=30-80 K for the 
estimates in Table O 

The physical parameters listed in Table [6] can be calculated from two independent 
datasets for each core (SMA 1.3 mm and CARMA 1.4 mm for Gil. 92-0.61, SMA 1.3 mm 
and CARMA 3.4 mm for G19.01— 0.03). For each compact millimeter continuum source 
in Table El the mass estimate derived from the lower resolution dataset is greater than 
that derived from the higher resolution dataset. Conversely, a larger beam- averaged column 
density and volume density are calculated from the higher resolution data. These trends are 
consistent with the lower-resolution data being more sensitive to emission on larger spatial 
scales. We note that the mass estimates derived from the dust continuum emission include 
only circum(proto) stellar material, and not the mass of any protostar or ZAMS star that 
has already formed within a compact core. 

For comparison. Table [5] presents estimates of N(H2), n(Il2), and Mgas derived from the 
best-fit source size and CH3CN column density for the hot cores G11.92— 0.61-MMl and 
G19.01-0.03-MM1. Estimates are presented for CH3CN/H2 abundances of 10"^, 10-^ and 
10~^. Values for the abundance of CH3CN/H2 in hot cores reported in the literature span 
at least an order of mag nitude, from ^lO^^-lO"'^ flRemiian et al.ll2004l : IZhang et al.l 120071 : 
Bisschop et al.l 120071 ). Lower abundances (~10^^) may also be possible even at relatively 
high temperatures (>100 K) in massive hot cores , depending on the warm-up timescale 
driving the gas-grain chemistry (iGarrod et al.ll2008l ). Given the uncertainty in the CH3CN 
abundance, the gas mass estimates derived from the CH3CN emission (Table [5]) and from 
the millimeter dust continuum emission (Table [6]) are broadly consistent. 



4.1.4- Nature of the Continuum Sources 

In summary, the millimeter continuum sources G11.92— 0.61-MMl, G11.92— 0.61-MM2, 
and G19.01— 0.03-MMl are dominated by thermal dust emission. The circum(proto) stellar 
gas masses of these cores range from ~8-62 M0 (based on the SMA data, resolution ~3x2"). 
Gil. 92— 0.61-MMl and G19.01— 0.03-MMl are hot cores, with derived gas temperatures of 
166± 20 K and 144± 15 K, respectively. SED modeling indicates that a central (proto)star of 
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~10 M0 is present within the G19.01— 0.03-MMl core. The properties of individual members 
of the Gl 1.92—0.61 (proto) cluster cannot be constrained by this method, as the sources 
MMl, MM2, and MM3 are unresolved in available data at wavelengths < 1.3 mm. However, 
the bolometric luminosities of G19.01— 0.03 and of the Gil. 92— 0.61 (proto) cluster as a whole 
are comparable (~10^ Lq). The nature of G11.92— 0.61-MM3 is less clear. In principle, an 
HCHII region undetected in previous observations could account for the majority of the 
G11.92— 0.61-MM3 mm flux density ( §4.1.3p . but additional observations at cm wavelengths 
are needed to investigate this possibility. If the mm flux density of G11.92— 0.61-MM3 is 
dominated by dust emission, the compact gas mass is ~3-9 M0, the lowest of the observed 
cores. The relative evolutionary states of the members of the Gil. 92— 0.61 (proto) cluster, 
and of Gil. 92-0.61 and G19.01-0.03, are discussed further in gXl 

Based on the SMA 1.3 mm data, the total mass in compact cores is ~37-94 Mq in 
Gil. 92— 0.61 and ~12-16 M© in G19.01— 0.03. Additional low-mass sources may also be 
present, but undetected in our observations; the 5cr sensitivity limit of the S MA data corre 



spends to a mass limit of a few Mq for moderate dust temperatures (Table [6]). ISchuUer et al. 



( I2OO9I ) calculate a mass for the larger-scale (40x34") G19.01-0.03 gas/dust clump of 1070 
Mq, based on ATLASGAL 870 fj,m data and an NH3 Tkin of 19.5 K. This suggests that only 
~1% of the total mass is in the compact core we observe with the SMA, and a considerable 
reservoir of material is in an extended envelope that is mostly resolved out in the continuum 
as in the ^^CO line emission ( §3.2.4p . The compact cores in the Gil. 92— 0.61 protocluster 
co nstitute a larger pe rcentage of the total mass reservoir. From the 850 /im SCUBA flux (12 
Jy, IWalsh et al.ll2003[ ). we estimate a total mass for the clump of ^780 M© for Trf„,,st =20 K 



(typical of the NH3 temperatures reported for ATLASGAL s ources bv ISchuUer et a 



and K850Atm=2.2 (interpolated from the values tabulated by lOssenkopf fc Henning 



20091) 



19M). 



Based on this estimate, the compact SMA cores in Gil. 92— 0.61 comprise ~5-12% of the total 
mass, with a remaining large-scale gas reservoir of several hundred M© for the Gil. 92— 0.61 
(proto) cluster. 

Single dish surveys of massive star forming regions have revealed spectrosco pic signatures 



of parsec-scale infall in cluster forming environments (e.g. iWu &: EvansI 120031 ). In addition, 
new high resolution observations of the G20.08— 0.14 N cluster detect infall at the scale 
of both cluster forming cl umps and massive star for ming cores, all part of a continuous, 
hierarchical accretion flow ( iGalvan-Madrid et al.ll2009l ). Recent simulations also indicate the 
importance of accretion from large-scale gas reservoirs in massive star and cluster forniation , 
particularly for d etermining the flnal stellar masses (jSmith et al.l l2009l : iPeters et al.l I2OIOI : 
Wang et al.ll2010l ). Since the presence of an active outflow indicates ongoing accretion, the 
masses of the members of the Gl 1.92—0.61 (proto) cluster may grow signiflcantly with time. 
For G19.01— 0.03, the SED modeling is consistent with a central YSO of mass ~10 M© that 
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is actively accreting at a rate of ~10~^ Mq year~^. This central source is associated with a 
compact gas and dust core of mass ~12-16 Mq. However, with a substantial (~1000 Mq) 
extended reservoir of material from which to draw, the final mass of G19.01— 0.03 may be 
substantially higher. 



4.2. Outflows 



A single dominant bipolar molecular outflow is associated with each of our targeted 
EGOs. These outflows are traced by high- velocity, well-coUimated ^^C0(2-l) and HC0'*"(1- 
0) emission. In both EGOs, the red and blue outflow lobes clearly trace back to a driving 
source identified with a compact mm continuum core (Figs. |6l [T3|) . This relative clarity is 
somewhat unusual. In many massive star-forming regions, multiple outfiows are observed, 



with c omplex kinematics th a t can make it difficu lt to identify driving source(s) (jZhang et al. 



20071 : IShepherd et al.l 120071 : iBrogan et al.l l2009l ). Indeed, si nce YSOs of all m asses drive 
bipolar molecular outfiows during the formation process (e.g. [Richer et al.ll2000l ). one would 
expect multiple outfiows in a protocluster such as Gil. 92— 0.61. 

A second outfiow may indeed be present in Gil. 92-0.61. Blueshifted ^^00(2-1), HC0+(1 
0),and SiO(2-l) emission are present NE of the mm continuum cores, and redshifted emission 
to the SW ( §3.1.4p . This is opposite the velocity gradient of the dominant outfiow, and this 
emission may trace a second outfiow. If so, the driving source is likely the continuum source 
MM3, which is approximately equidistant between the two lobes (Fig.[6t; the possible second 
outfiow is most prominent at moderate velocities, see also §3.1.4p . Alternatively, the observed 
morphology may be attributable to orientation effects. An outflow near ly in the plane of 



the s ky may appear to have overlapping red and blue-shifted lobes (e.g. ICabrit fc Bertout 



19901 ). Another possible explanation is outflow precession. For an outflow axis close to the 
plane of the sky, precession can produce the appearance of i nversions between blue/red- 



shifted emission along the outflow axis (e.g. iBeuther et al.ll2008l ). such as the pattern seen in 



Gil. 92— 0.61. In addition, the ^^CO and HCO"^ data hint at the possible presence of a third, 
low-velocity outflow along a SE-NW axis. As shown in Figures [THHl moderately redshifted 
gas is present SE of the continuum sources, and moderately blueshifted gas to the NW (26.7, 
39.9, and 46.5 km s~^ panels). The interpretation of this emission as an outflow is, however, 
very uncertain. The moderate-velocity ^^CO emission appears to correlate with extended 
4.5 /im emission and 44 GHz CH3OH masers, but the HCO"'" emission (which is subject to 
less spatial flltering) is much more extended, suggesting confusion with the ambient cloud, 
and the SiO(2-l) emission (Fig. [9]) does not show the same velocity pattern. There is no 
clear evidence in our data for an outflow driven by the continuum source MM2. 
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4-2.1. Outflow Properties 

We estimate the physical properties of the molecular outflows in Gil. 92— 0.61 and 
G19.01-0.03 independently from the SMA i2C0(2-l) and the CARMA HCO+(1-0) data. 
As discussed in §3.1.41 and §3.2.41 outflow gas is confused with diffuse emission from the 
surrounding cloud at velocities near the source vlsr- This problem is particularly acute for 
^^CO, because of its high abundance. To avoid including contributions from the ambient 
cloud, we consider only high velocity gas in our estimates of the outflow physical properties 
(Table [7]). To further isolate the outflow gas, a polygonal mask is deflned for each red or 
blueshifted outflow lobe in Figures E] and [T31 The polygonal masks are drawn to encompass 
all outflow emission in the integrated intensity images of the high-velocity gas, and checked 
against the datacubes. The appropriate mask is applied to each channel in which the outflow 
dominates over emission from the ambient cloud. Assuming optically thin emission, the gas 
mass of the outflow is then calculated from 



M, 



out 



1.186 X 10-^Q{Tecc)e^^D^ f S^dv 



(4) 



where Mout is the outflow gas mass in Mq, T^x is the excitation temperature of the transition 
in K, Q(Te^.) is the partition function, Eupper is the upper energy level of the transition in K, z/ 
is the frequency of the transition in GHz, x is the abundance of the observed molec ule relative 



to H2, D is the distance to the source in kpc, and is the line flux in Jy. Following iQiu et al. 



(l2009l ). for ^^CO we adopt an abundance (x) relative to H2 of 10 an excitation temperature 



of 30 K, and a mean gas atomic weight of 1.36 (included in the constant in equation (1)). 
For HCO"*", we ado pt the same excitation temperature (Tp .T ,=30 K), and an abundan ce of 
10"^ relative to H2 JVogel et al.]ll984l : iRawhngs etaP 120041 : Ixiaassen fc Wilson! 120071 ) . We 
use Q(30 K) = 11.19 for ^^CO and Q(30 K) = 14.36 for HCO+, interpolati ng from the values 



provi ded in the Cologne Database for Molecular Spectroscopy (CDMS. iMiiller et al.l 12001 
2OO5I ) and /i2S=0.02423 debye^ for ^200(2-1) and /x^S=15.21022 de bye^ for HCO+(1-0) from 



the Splatalogu^ spectral line database. Following Qiu et al. (2009), we estimate the outflow 
momentum and energy using 

Pout = ^Mout{^v)Av (5) 

and _^ 

Eout = -^Mout{Av){Av)^ (6) 



where MoMf(Av) is the outflow mass in a given channel and Av =\vcenter,channei—VLSR\- ^oi 
these calculations, we adopt fL5_R=35 km s~^ for Gil. 92— 0.61 and fL5i?=60 km s~^ for 



%ttp: / /www. splatalogue.net / 
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G19.01— 0.03. We estimate the dynamical timescale from tdyn = Loutfiow/vmax where the 
length Loutfiow and the maximum velocity Vmax are determined separately for the red and 
blue lobes of each outflow (Table [7]). In estimating Loutfiow, we measured the extent of 
the red/blueshifted emission from the driving mm continuum source. For Gil. 92— 0.61, we 
assumed that the main outflow was driven by MMl, and the possible second ("northern") 
outflow by MM3. Using the dynamical timescales, we also estimate the mass and momen- 
tum outflow rates, Mont = Mout/tdyn and Pout = Pout/tdyn- FoT each outflow, the outflow 
parameters are listed in Table [TJ along with the velocity ranges used. For G19.01— 0.03, the 
"NE blue clump" (Table [7]) is the easternmost knot of blueshifted ^^CO emission (Fig. [T3k). 
This knot is offset from the main ^^CO jet, and a separate mask was defined for it. However, 
the HCO"'" and SiO morphology indicate that this ^^CO emission is likely part of the outflow, 
so we include it in our estimates of the outflow properties. 

Several salient points are reflected in Table [3 (1) channels nearest the systemic velocity 
disproportionately affect the outflow mass estimates; (2) the estimates derived from HCO"*" 
and ^^CO observations differ by approximately an order of magnitude; and (3) there is 
considerable uncertainty in the estimate of the dynamical timescale, and hence of the mass 
and momentum outflow rates. Each of these points is discussed in more detail below. 

As noted above, dominant, unconfused outflow emission is the primary criterion for 
choosing the velocity (channel) ranges over which to integrate outflow mass, momentum, 
and energy. In general, column densities are highest near the systemic velocity of a cloud. 
As a result, estimates of outflow mass and other properties are extremely sensitive to 
how closely the velocities considered approach the vlsr, e.g. to the minimum value of 
Ay=\y center,channei—VLSR\- For thls reasou, where practicable we choose velocity ranges such 
that min.(Avf,i„e)~min.(Avred) and used the same velocity range for ^^CO and HCO+ mass 
estimates. In G19.01— 0.03, it is possible to follow the outflow much closer to the systemic 
velocity in HCO+(1-0) than in ^^C0(2-l), with minimal confusion from ambient diffuse gas. 
As an illustrative example, in Table [7] we present estimates of the G19.01— 0.03 outflow 
properties derived from HCO^ using velocity ranges beginning ~±6km s~^, ±8km s~^, and 
±10km s~^ from the systemic velocity. The difference in the estimated total outflow mass 
(and consequently in Mout) is about a factor of 2. The estimates of the outflow momentum, 
energy, and momentum outflow rate are less severely affected because the channels in ques- 
tion are near the systemic velocity (low Av). In Gil. 92— 0.61, the situation is complicated 
by the possible second outflow, so we restrict the HCO+ mass estimates to the same velocity 
range considered for ^^CO. 

For a given outflow lobe, the mass estimate derived from IICO^(l-O) is roughly an order 
of magnitude larger than that derived from ^^C0(2-l). There are two possible explanations 
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for this discrepancy: spatial filtering a nd uncerta i n HC O"*" abundance. For the massive 
outflow in G240.31+0.07 (D=6.4 kpc), lOiu et al.l toO^ found that their SMA compact 
configuration ^^CO observations recovered <10% of the single-dish flux at line center, ~80% 
in the "line wing" (Av~13 km s~^), and nearly 100% at more extreme velocities (Av> 15 
km s~^). Our CARMA HCO"'" data are m ore sensitiv e to la rger-scale emission than our SMA 
^^CO data, and the linear resolution of the lQiu et al.l (120091 ) SMA observations is comparable 
to that of our CARMA EGO data. It is plausible that our CARMA observations are picking 
up outflow emission on larger spatial scales, to which the SMA is insensitive. In this case, 
the outflow parameters estimated from the HCO^ emission would be more representative 
of the true outflow properties. If, however, the IICO+ abundance in our target EGOs is 
enhanced above our adopted value of 10~^, this could lower the mass esti mates from the 
HCO" *" emission into better agreement with those from ^^CO. Modeling by iRawlings et al. 
( I2OO4J ) found best-fit HCO"*" abundances of 10~^° in the envelope, 10~^ in the jet/cavity, 
and 10~^ in the boundary layer, though the models were optimized for the low-mass source 
L1527. Moderate optical depth corrections would also increase the ^^CO mass estimates (see 
below), bringing them into better agreement with those from HCO^. 

There is considerable ambiguity in the determination of the dynamical timescale, par- 
ticularly for asymmetric and/or clumpy outflows such as those observed towards our target 
"""" ' " " ' ' - ■ ^ " 1 interferometric outflow studies (e.g. 

Lout flow / 'Vr, 



Qiu & Zhang 


2009; 


Qiu et al. 


2009) 



and calculate tdyn inde- 



pendently for the red and blue lobes of each outflow. As Table [7] illustrates, these estimates 
can differ significantly. Some single-dish studies instead use tdyn = R/V ^ where R is the dis- 
tance between the peaks of the red and blue outflow lobes, and V is the mean out flow velocity, 



calcu lated as the outflow momentum divided by the outflow mass (P/M, e.g. IZhang et al. 



2OO5I ). Both approaches assume, however, that an outflow (or an outflow lobe) can be well- 



characterized by a single velocity. The G19.01— 0.03 outflow is clumpy, and characterized by 
discrete knots of high velocity gas. If tdyn = L knot I v max, knot is calculated individually for each 
of the three blueshifted knots in Figure [T3k. the values are (from north to south) ~3600 years, 
1100 years, and 600 years. This range suggests the limitations of attempting to evaluate the 
age of a flow using a single velocity. Estimating idyn is further complicated by the unknown 
effects of ambient density and inclination angle. The expressions for idyn above assume free 
expansion. Finally, the inclination of the outflows is unknown. The extended morphologies 
of the 4.5 iim. emission and the high- velocity molecular gas in our target EGOs suggest that 
the outflows may lie near the plane of the sky (Figs.[6| [T3|) . However, intermediate inclination 
angles are also plausible, since th e red and blueshifted o utflow lobes are spatially distinct, 
particularly in G19.01— 0.03 (e.g. ICabrit fc Bertoutlll990l ). and very high velocity gas (~60 
km s~^ from vlsr in Gil. 92— 0.61 and >100 km s~^ from vlsr in G19.01— 0.03) is observed. 
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Table [7] presents estimates of outflow parameters without correction for inclination, and for 
^=10°, 30°, and 60°, where 9 is the inclination angle of the outflow to the plane of the sky. 
In the extreme case of ^=10°, correcting for inclination increases the estimated yiout and 
Poui by a factor of ~6, and Vout and 'Ejout by a factor of ~30. For an intermediate inclination 
6^=30°, the increases are more moderate: a factor of ~2 for yiout and Pout, and ~4 for Vout 
and 'Eiout- For outflows in which the red and blue lobes give very different estimates of tdyn, 
Table [7] presents estimates of yiout and Pout for the outflow as a whole using an intermedi- 
ate timescale. In calculating Viout and Pout from the HCO^ data we adopt the dynamical 
timescales calculated from ^^CO, since the highest velocity outflow gas extends beyond the 
limited velocity coverage of the CARMA observations. 

In general, our estimates of outflow mass are lower limits, and likely extreme lower 
limits. As a result, the other physical parameters (which depend on the outflow mass) will 
also be underestimated. There are three main contributing factors: (1) extended emission 
missed by the interferometers; (2) outflow emission near the systemic velocity excluded by our 
conservative velocity cuts; and (3) the assumption of optically thin emission. Our estimates 
of the outflow mass assume optically thin emission in both ^^C0(2-l) and HCO^(l-O). While 
this assumption is likely valid for HCO"'", it is more problematic for ^^CO, and some recent 
interferometric outflow studies have made detailed corrections for ^^CO optical depth (e.g. 



Qiu et al.ll2009l ). Because we were conservative in selecting the velocity ranges over which 
to calculate outflow parameters, significant (>5cr) ^'^00(2-1) emission is detected in only 
one channel that contributes to the estimates presented in Table [3, for one outflow: the 
main outflow in Gll.92-0.61 (v=48.1 km s'^). If we calculate the ^^CO/^^CO line ratio 
at the ^^CO peak in this channel, the implied optical depth correction fac tor is ^6.7 for a 
^^CO /^'^CO abundance ratio of 40 (for a Galactocentric distance of 4.7 kpc, IWilson fc Rood 



19941 ). Applying this factor would increase the contribution of this single channel to the 
outflow mass from ~0.07 Mq to ~0.5 M©, the mass of the red outflow lobe from ~0.2 to 
~0.6 Mq, and the total mass of the outflow from ~0.8 M0to ~1.2 Mq. Applying this single 
correction factor, however, would likely result in an overestimate. The correction factors 
derived at two other p ositions in the out f low w ith detected ^^CO emission are more modest 
(~3 and 4.5, see also ICabrit fc BertoutI Il990l ). The signal-to-noise of the ^'^CO data are 



not sufficient to accomodate attempting an opacity correction as a function of position, 
particularly given the overwhelming lack of detected ^'^CO emission in the other channels 
considered. By assuming optically thin emission, our mass estimates based on ^^CO will 
definitively be lower limits, without the ambiguity of possibly overcorrecting. 

We do not estimate outfiow parameters from the SiO(2-l) emission observed with 
CARMA because of the uncertainty in the SiO abundance in the emitting region. Values in 
the literature for the SiO abundance in molecular outfiows and massive star-forming regions 
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vary over at least two orders of magnitude (from ~10 ^ to 10 ®, iPineau des Forets et al. 



19971 : ICaselli et al.lll997l : ISchilke et al.lll997| ). Models indicate that the SiO abundance de- 



pends sensitively on the shoc k conditions (including velocity, ambient density , and time 
since the passage of the shock, iPineau des Forets et al.lll997l : ISchilke et al.l 119971 ) . which are 
not constrained by our single-transition SiO observations. Since our CARMA data show 
that the SiO emission is extended well beyond the beamsize of the JCMT SiO (5-4) spectra 
( §3.1.H§3.2.ip . we cannot constrain the physical conditions in the SiO emitting gas. 



4-2.2. Comparison with Other Objects 



Properties of outflows from MYSOs reported in the literature, based on high angular 
resolution observations, range over several orders of magnitude. As indicated by the discus- 
sion above, some of this range may be attributable to differences in spatial filtering and to 
the (large) uncertainties inherent in assuming tracer abundances and corr ecting (or not) for 
optical depth and inclination effects. At the low end, lOiu fc Zhangl (120091 ) calculate Mout of 
0.22 Mq, Pout of 4.9 Mq km s-\ Mont of 10"^ Mq jr-\ and Pout of 2.2x10^^ km s'^ 
yr~^ based on SMA ^^CO data for the extremely high velocity outflow in HH80-81 (D=1.7 



kpc) . At the hig 



kpc: IZhaiig et al 



Oiu et al 



2009 



1 end , outflow masses of several tens (27 Mq, IRAS 18566+0408, D=6.7 
20071) to > 100 M^ (98 Mq, G240.31+0.07, D=6.4 kpc; 124 Mq, Orion-KL: 



Beuther fc NissenI |2008| ) have been reported. These st udies, however, use 



tracers with uncertain abunda nce in outflows f SiO: IZhang et al.l 



dish and interferometric data (jOiu et al.l l2009l : iBeuther fc Nissen 



2003), or combine single 
2008i ). Also, except for 



Orion-KL, the estimated dynamical timescales for these more massive outflows are longer 
(>10^ years), so the estimated mass outflow rates, Mout, are still only a few x 10^^ Mq 
yr^^. The estimated parameters of the molecular outflows in our target EGOs (Table [7]) 
are roughly in the middle of the range reported in the literature. The main outflow in 
Gil. 92— 0.61 and the outflow in G19.01— 0.03 have broadly similar characteristics: each has 
tdyn of a few xlO^ years, and (based on the HCO"*" data) Mout of a few Mq, Pout of a few 
hundred Mq km s~^, Eout of tens to a hundred xlO"^^ ergs, Mo«t of a few xlO~'^ Mq yr~^, 
and Pout of a few hundredths to one Mq km s^^ yr~^ (the estimates of Eout and Pout are 
most severely affected by the uncertainty in the inclination angle). These parameters are 
generally comparable to those for the total high-velocity gas (attributed to three separate 
outflows) in the massive star -forming regio n IRA S 17233-3606 derived from high-resolution 
SMA ^^CO observations by iLeurini et al.l (120091 ) (opacity correction applied), though for 
IRAS 17233-3606 the estimated dynamical timescale is somewhat shorter (~300-1600 years) 
than for our target EGOs. The EGO outflow properties are also quite similar to those of 
the outflows in the AFGL5142 protocluster, as estimated from OVRO HCO+(1-0) and SMA 
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C0 (2-l) observations (particularly accoun ting; for the different assumed HCO^ abundance, 



IQ- ^: iHunter et al.lll999l : IZhang et all 120071 ) . As noted in g3X3l the frequency coverage of 
the IZhang et al.l (120071 ) SMA data is comparable to that of our observations, and the SMA 
sp ectrum of AFGL5142 MM2-the probable driving source of the north-south outflow studied 
by iHunter et al.l (jl999l )-is very similar to that of G11.92— 0.61-MMl. Even the least mas- 
sive outflow observed towards our target EGOs (the "northern outflow" in Gil. 92— 0.61) 
has values of Mout, ^out, etc. at least an order of magnit ude greater than thos e typical of 
low- mass outflows observed at high angular resolution(e.g. Arce Sz Sargentl[2006l ). 



Several large scale single-dish surveys of molecular outflows have shown correlations 
between the properties of the outflow an d those of the dr i ving source (in particular its 
bolometric luminosity and core mass, e.g. ICabrit fc BertoutI Il992l : IShepherd &: Churchwell 
19961 : lHunterlll997l : iBeuther et al.ll2002l ) , though other rec ent studies have found considerable 



scatt er and weak or no evidence of any correlations (e.g. iRidge fc Moorell200ll : IZhang et al. 



20051 ) . The applicability of these relations to parameters derived from interferometric obser- 
vations is unclear, since, as discussed above, interferometers resolve out a significant fraction 
of the extended emission and so underestimate the outflow mass and other parameters. The 
mass and momentum outflow rates for our target EGOs (Table [7]) do generally agree, within 
the considerable scatter, with the Mout and Pout expected for a driving source of luminosity 
~10^ Lq {Mout ^a few x 10"^ to a few x 10"^ M^r, year-\ P„„i^a few x 10"=^ to lO"^ M,^ 
km s~^ year"^: ICabrit fc BertoutI Il992l : IShepherd fc GhurchwelM 1 19961 : iBeuther et all 12002 



Zhang et all 120051 ). For Gil. 92-0.61 and G19.01-0.03, the mass in the outflow (as derived 
from the HCO^ data. Table [7]) is roughly comparable to the circum(proto) stellar gas mass 
of the driving compact millimeter core (Table |6]), and, for G19.01— 0.03, to the mass of the 
central (proto)star inferred from the SED modeling ( §4.1.ip . It has been suggested for some 
time that the ma ss in MYSO outflows is accumulated from the larger-scale environment (e.g. 
Churchwelll Il997l ) . Notably, the properties of the outflows in our target EGOs are consistent 
with the single-dish relations with respect to "core mass" only if the total masses of the 
single-dish clump (as opposed to the masses of the compact cores resolved with the SMA 
and CARMA) are considered. 



4.3. Diversity within the EGO Sample 

Our high resolution mm observations suggest considerable diversity within the EGO 
sample in the clustering properties and evolutionary states of the outflow driving sources. 
Since extended 4.5 fim emission specifically targets a population with ongoing outflow activ- 
ity and active, rapid accretion, the range in source properties is of interest in the context of 
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recent theoretical work on feedback effects in cluster-scale models of massive star formation 
(e.g. 



Krumholz fc Matzned l2009l : lBatell2009l : IWang et al.l 120101 : iPeters et al.ll2010f ). 



4.3.1. Multiplicity 

The EGO Gl 1.92—0.61 was chosen as a target for high resolution mm observations in 
part because its 24 fim morphology and two associated 6.7 GHz Class 11 CH3OH maser 
spots indicated the possible presence of multiple MYSOs. Our high-resolution SMA and 
CARMA data indeed reveal three compact millimeter continuum sources associated with 
the EGO. The clustering scale of these cores is ~0.1 pc, comparable to that in S255N, 



a ma ssive protocluster that is also associated with a 4.5 /im outflow ( ICyganowski et al. 



20071 ). No additional structure is seen in the CARMA 1.4 mm image (resolution 1'.'44 x 0'.'87 
~0.03x 0.02 pc ~5500 x 3300 AU) as compared to the lower resolution SMA 1.3 mm image 
(resolution 3'.'2 xl'.'S ~0.06 x 0.03 pc ~12000 x 6800 AU). Further structure may, however, 
be suggested by the complicated kinematics of the ^^CO, HCO+, and particularly the SiO 
emission. In NGC6334I(N), a tight (~800 AU separation) binary is thought to cause the 
precession of the molecular outfl ow: both coniponen ts of the binary are contained within a 
hot core of diameter ~1400 AU (IBrogan et al.ll2009l ). If outflow precession is responsible for 
the complex velocity structure in Gil. 92— 0.61, it would suggest that the outflow driving 
source, G11.92— 0.61-MMl, could be an unresolved (proto)binary. 

In contrast to Gil. 92— 0.61, our interferometric observations of G19.01— 0.03 reveal only 
a single continuum source. For G19.01— 0.03, our SMA 1.3 mm continuum image provides the 
highest resolution: 3'.'2xl'.'7 ~0.07x0.03 pc ~13400 x7100 AU. This is sufficient to resolve 
clustering at the scale seen in Gil. 92— 0.6 1 and S255N, though not proto- Trapezia (such 
as NGC63341 and IfNV lHunter et al.l 120061 ). The CH3CN(12-11) spectrum of G19.01-0.03- 
MMl is unusual, in that the k=0, k=l, and k=2 components are nearly equal in strength, 
and the k=3 component is even stronger (Fig. |2T|) . This is likely indicative of either un- 
resolved multiplicity or unresolved density/temperature gradients, but the current data are 
insufficient to distinguish between these scenarios. 

If G19.01— 0.03 is truly a single (proto)star, it would be an unusual example of a massive 
star forming in isolation. Some recent models suggest that an appa rently lone MYSO could 
also be indicative of the very earliest stages of cluster formation. iPeters et al.l ( I2OIOI ) find 
that the most massive star begins to grow early in the cluster formation process, and fin- 
ishes accreting while its surrounding cluster is still forming. Accretion is clearly ongoing in 
G19.01— 0.03, and our mm observations provide several other indicators of its youth ( §4.3.2p . 
Higher resolution and more sensitive (sub)mm observations of G19.01— 0.03 are needed to 
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determine whether sub- 10000 AU clustering or low-mass (proto) stars are present. 



4-3.2. Relative Evolutionary State 

The three millimeter continuum sources in the Gil. 92— 0.61 (proto) cluster exhibit a 
range of MIR-mm properties and maser associations, suggestive of a range in mass and/or 
evolutionary state. The strongest mm continuum source, G11.92— 0.61-MMl, is a hot 
core, and the driving source of the dominant bipolar outflow in the region. In contrast, 
G11.92— 0.61-MM2 is devoid of millimeter-wavelength line emission, and shows no evidence 
for associated outflow or maser activity. The gas mass of G11.92— 0.61-MM2, calculated from 
its arcsecond-scale thermal dust emission, is ~17-42 M©, sufficient to form an intermediate to 
high-mass star. The compact gas mass associated with G11.92— 0.61-MMl is smaller, ~5-13 
Mq. This does not, however, include the mass of central (proto)star(s), which the associated 
MIR and maser emission indicate are likely present. MM2, which exhibits no signs of maser 
or outflow activity, may be a massive protostellar core. High-resolution observations in a 
thermometer of cold, dense gas (e.g. NH3) would help to determine whether this source is 
internally heated by a central YSO or externally heated by feedback processes in the cluster 
environment. 

The evolutionary state of G11.92— 0.61-MM3, the weakest millimeter continuum source 
in the region, is somewhat ambiguous. It is associated with 6.7 GHz Class II CH3OH maser 
and 24 /xm emission, and with the strongest 8.0 fim emission in the region (Fig. [T^,b). 
The nature of the 8 fim counterpart is unclear, as it appears somewhat extended and it is 
possible that H2 line emission contributes significantly to the broadband fiux. Very little 
compact mm molecular line emission is detected coincident with MM3, and no emission in 
high excitation (Eupper>100 K) molecular lines. Like MMl, the compact gas mass associated 
with MM3 is relatively small (~2-6 Mq), but t he MIR and mas er emission suggest that 



a central (proto) star has likely already formed. iDe Buizerl ( 120061 ) suggested that Class II 
CH3OH maser emission may be excited in the MIR-bright walls of an outfiow cavity. The 
morphology of the 24 /xm emission, however, appears inconsistent with MM3 being a hotspot 
on an outfiow cavity wall. Another possibility is that MM3 is more evolved, and so has 
dispersed most of the molecular gas in its immediate vicinity. This would be consistent with 
its MIR properties. We also note that MM3 is the only compact millimeter continuum source 
for which the arcsecond-scale 1.4 mm emissio n observed with CARMA could conceivably be 



due to a HC HII region that fell below the ICyganowski et al. J2009h 44 GHz continuum 



detection limit ( §4.1.3p . High resolution observations over a range of cm wavelengths are 
required to constrain the evolutionary state of MM3. 
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Several lines of evidence suggest that G19.01— 0.03 is younger than the Gil. 92— 0.61 
massive star-forming region. Compared to G11.92— 0.61-MMl, the millimeter molecular 
line emission from G19.01— 0.03-MMl is sparse and weak. The relative lack of chemical 
complexity observed towards G19.01— 0.03-MMl, and its co oler derived temperature, ar e 
both consistent with a less-advanced hot core chemisty (e.g. Ivan Dishoeck fc Blakd Il998l ). 
There have been a number of suggestions in the literature that the relative abundances 
of various sulfur-bearing species may provid e chemical clocks for hot cores (e.g. ICharnley 



19971 : iHatchell et al.lll998l : iHerpin et al.ll2009l ). Transitions of several species invoked in these 
models fall within the frequency coverage of our SMA observations, namely SO, OCS, and 
SO3. Unfortun ately, the most diagnostic abundance ratios are those relative to H2S (e.g. 
Charnleylll997l ). which does not have transitions within our SMA bandwidth. 



The calculation and comparison of abundances for the sulfur-bearing species in the 
EGO hot cores must be regarded with considerable caution. Such comparisons rely on an 
assumption that all relevant emission arises from the same excitation conditions in the same 
physical volume of gas. From our data, it is evident that in Gil. 92— 0.61 the outflow con- 
tributes to the excitation of SO emission, and even the compact OCS emission may have 
an outflow contribution ( §3.1.3p . Considering the two-component temperature model re- 
quired to fit the CH3CN emission, it is also plausible that the low-excitation SO emission 
arises from a larger volume than the higher-excitation OCS emission (both unresolved by the 
SMA beam). Similarly, the thermal dust emission (from which we infer N(H2)) may arise 
from a different volume than the molecular line emission. If we nonetheless compute the 
OCS and SO abundances for G11.92-0.61-MM1 assuming Te^=166 K and N(H2)=2. 7x10^3 
cm~^ (beam-averaged column density for Tdust=^QQ K and the SMA beam), we find abun- 
dances of ~1.9 X 10-^ and 1.3 x 10"^ relative to H2, respectively. For G19.01-0.03-MM1, 
(Te^=115 K, N(H2)=3. 5x10^3 cm'^. Table [6]), we calculate abundances of ^3x10^^ for 
OCS and ~6xl0'^° for SO. These estimates correspond to a [SO/OCS] ratio of ~0.2 for 
G19.01-0.03 and ~0.7 for Gil 92- 0.61-MMl, consistent with G19.01-0.03 being younger 
and less evolved (IHerpin et al.ll2009l ). Other, more global indicators also point to the rela- 
tive youth of G19.01— 0.03 compared to Gil. 92— 0.61. Since MYSOs dissipate their natal 
envelopes as they grow, the presence of an extended envelope is suggestive of youth. Compact 
mm core(s) account for a smaller fraction of the total clump mass in G19.01— 0.03 than in 
Gl 1.92-0.61 (SX3]), consistent with indications from the ^^CO emission and SED modeling 
that G19. 01— 0.03-MMl is still embedded in a massive (~1000 Mq) large-scale envelope. 
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4.3.3. Future Work 

These initial results demonstrate the potential of the EGO sample for probing the im- 
portance of protostellar feedback in the formation of massive stars and star clusters. Recent 
theoretical work has just begu n to include realist ic feedback effects in cluster-scal e models, in- 



cluding protostellar outflows (IWang et al.ll2010l ). photoioniz ation/HII r egions ( iPeters et al. 
20 id : iKrumholz fc Matznerl [20091 ) . and radiative feedback jBatelboOoh . However, current 



models do not include all feedback mechanisms, and hence do not address the question 
of which mechanism(s) are most important at which stages of (proto) cluster formation, or 
of how these mechanisms interact. EGOs are an outflow-selected sample. Hence, charac- 
terizing their outflows as well as other possible forms of feedback (ionized jets/winds, gas 
heating) will help to establish whether there is an outflow-dominated stage of (proto) cluster 
formation. Addressing this question in a statistically meaningful way will require high res- 
olution (sub)mm wavelength data similar to those presented here for a wide range of EGO 
sources. High-resolution cm continuum data are also necessary to constrain the presence and 
physical properties of ionized gas, while cm-wavelength line observations (e.g. of NH3) are 
needed to constrain gas temperatures and densities in cool cores that lack mm-wavelength 
line emission. Accumulating uniform datasets for large samples is essential for comparisons 
of different objects. The capabilities of the EVLA and ALMA put such surveys within reach, 
and EGOs will be promising targets for these instruments. One fortunate characteristic of 
EGOs is the probable 229.759 GHz CH3OH Class I maser line which we detect toward both 
objects. These are, to our knowledge, the flfth and sixth reported examples of this maser in 
the literature, although higher resolution observations are still required to unambiquously 
conflrm the maser nature of the emission based on the brightness temperature. As with the 
44 GHz Class I masers, these features may be common in massive star-forming regions. If so, 
they will be of interest as self-calibration targets to help enable future very high-resolution 
1.3 mm observations (e.g. with ALMA). 



5. Conclusions 

Our high-resolution millimeter observations of two EGOs unambiguously show that they 
are young MYSOs driving massive bipolar outflows. The spatial coincidence of high velocity 
^^C0(2-l) and HCO^(l-O) emission with the extended 4.5 /im lobes supports the outflow 
hypothesis for the 4.5 fim emission. A single dominant outflow is identifled in each EGO, 
with tentative evidence for multiple outflows in one source (Gil. 92— 0.61). The morphology 
and kinematics of the SiO(2-l) emission differ from the other outflow tracers in that some 
of the strongest red and blueshifted features are offset from the extended 4.5 /im emission. 
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and may trace the impact of outflow shocks on dense gas in the surrounding cloud. The 
morphology of the high- velocity gas with respect to 44 GHz Class I CH3OH maser emission 
further solidifies the association of this type of maser with outflows. Anomalously intense 
and narrow components of 229.759 GHz CH3OH emission are also detected in the outflow 
lobes from both objects, suggesting additional Class I maser activity. The outflow driving 
sources appear as compact cores of millimeter continuum emission and dense gas, including 
the hot core molecules CH3OH, CH3CN and OCS. Coincident with 22 GHz water maser 
emisson, G11.92— 0.61-MMl shows considerably richer and stronger hot core line emission 
than G19.01— 0.03-MMl, consistent with its warmer temperature derived from the multi- 
transition analysis of the CH3CN and CH3OH emission (166±20 v. 114±15 K). Both hot 
cores exhibit 24 /xm and 70 fim emission in MIPSGAL images and contain 6.7 GHz Class II 
CH3OH masers, all consistent with their identification as MYSOs. 

Our observations also reveal considerable diversity within the EGO sample. Although 
observed at the same spatial resolution, G19.01— 0.03 appears as a single MYSO while 
Gil. 92— 0.61 resolves into a cluster of three compact dust cores. In addition to the differ- 
ence in multiplicity, several other factors point to G19.01— 0.03 being in a earlier evolutionary 
stage: SED modeling, the relative weakness of its hot core emission, and the dominance of 
the extended envelope of molecular gas. In contrast, Gil. 92— 0.61 appears to have already 
formed a protocluster whose members span a range of ages - one is a hot core and two are 
almost entirely devoid of line emission. These initial results demonstrate the potential of 
the EGO sample for probing the importance of protostellar feedback in the formation of 
massive stars and star clusters. The future capabilities of the EVLA and ALMA will enable 
uniform surveys of a statistically meaningful number of regions which will enable the relative 
importance of outflows, photoionization, and radiative feedback to be assessed. 
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Table 1. Observational Parameters 



Obs. Date Syn. Beam Prim. Beam Cont. 



Line Observations 



FWHP 



mjy beam ^ km s ^ mjy beam ^ 



Calibrators 
Gain Bandpass Flux 



Gll.92-0.61 



1.3 mm SMA 

2008 June 23 3.2x1.8 55 

1.4 mm CARMA 

2008 April 25 1.44x0.87 31 

3.4 mm CARMA 
2008 July 29 6.8x5.1 80 



3.5 



4.3 



1.1 



1.7 



100 J1733-130,J1911-201 30454.3 Uranus 



J1733-130,J1911-201 30454.3 J1733-130'' 



10 J1733-130,J1911-201 30454.3 Uranus 



G19.01-0.03 



1.3 mm SMA 

2008 June 23 3.2x1.7 55 

3.4 mm CARMA 

2008 July 30 5.7x5.1 80 



3.5 



0.5 



1.1 



1.7 



100 J1733-130,J1911-201 30454.3 Uranus 



10 J1733-130,J1911-201 30454.3 Uranus 



^Hanning smoothed 
''Assuming 8(1.4 mm)= 2.69 Jy 
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Table 2. Observed Properties of Millimeter Continuum Sources 





J2000.0 Coordinates'' 


Peak Intensity*^ 


Integ. Flux Density'' 


Size'^ 


Source 


, „ (h m S-) 


5 {° ' ") 


mJy beam~^ 


mJy 


"x " [P.A. °] 




Gil. 92-0.61 


1.3 mm 


SMA<^ 










MMl 


18 13 58.103 


-18 54 20.19 


219 (6) 


338 (14) 


1.7x1.7 [126] 


MM2 


18 13 57.860 


-18 54 14.04 


136 (6) 


203 (14) 


1.9x0.8 [161] 


MM3 


18 13 58.21 


-18 54 16.2 


49 (6) 






Total 








590 (34) 




1.4 mm 


CARMA° 










MMl 


18 13 58.116 


-18 54 20.14 


96 (5) 


169 (13) 


1.26x0.76 [26] 


MM2 


18 13 57.860 


-18 54 14.17 


78 (6) 


121 (13) 


1.20x0.55 [36] 


MM3 


18 13 58.159 


-18 54 16.0 


25 (6) 


36 (13) 


1.59X <0.57 [40] 


Total 








326 (39) 




G19.01-0.03 


1.3 mm 


SMA* 










MMl 


18 25 44.790 


-12 22 45.86 


207 (3) 


275 (7) 


1.9x0.9 [78] 


3.4 mm 


CARMAs 










MMl 


18 25 44.80 


-12 22 46.8 


9.4 (0.9) 


27 (3) 


8.5x5.9 [134] 



''Centroid position determined by fitting a single two-dimensional Gaussian component to each source. 
The statistical uncertainty is indicated by the number of significant figures; absolute positional uncer- 
tainties are at least an order of magnitude larger (see [[2]l. 

''Estimated statistical uncertainties are given in parentheses. Unresolved sources are indicated by 
in the integrated flux density column. 

'^Deconvolved source size determined by fitting a single Gaussian component to each source. Unre- 
solved sources are indicated by 

"^Synthesized beam 3'.'2xl'.'8 (P.A.=59°) 

''Synthesized beam iy44x0'.'87 (P.A.=25'') 

•"Synthesized beam 3^2x1^7 (P.A. =63'') 

sSynthesized beam 5^7x5^1 (P.A. =-27°) 



Table 3. 1.3 mm Spectral Lines: Gil. 92— 0.61 



Fitted Line Parameters 



Species 


Transition 


Frequency 


Euppcr 


/k Intensity^ 


\T ^ 

V center 


Width'' 


/ S dv'' 






(GHz) 


(K) 


Jy beam~^ 


km s^-*- 


km s^^ 


Jy beam~^km 




MMl 


HCqNfV7 = lf~) 


24-23 


219.174 


452 


0.53 (0.05) 


33.6 (0.6) 


13.1 (1.3) 


7 4 d 0") 


C2H5CN 


24o oo-2.)0 01 


219.506 


136 


0.42 (0.05) 


36.6 (0.6) 


10.9 (1.4) 


4.8 (0.8) 




2-1 


219.560 


16 


1.6 (0.1) 


34.9 (0.1) 


4.3 (0.3) 


7.2 (0.7) 


HNCO 


IO9 Q-99 


219.734 


228 










HNCO 


IO9 S-99 7"^ 


219.737 


228 










HNCO 


lOn in-Qn q*^ 


219.798 


58 


1 10 ("0 06) 


35 2 (n 3) 


q q fO 6") 


11 5 Co 


H13CO 


3i,2-2i,id 


219.909 


33 


0.59 (0.03) 


35.1 (0.2) 


9.1 (0.5) 


5.7 (0.5) 


SO 


65-54'' 


219.949 


35 


2.62 (0.05) 


35.3 (0.1) 


10.5 (0.2) 


29.2 (0.9) 


CH3OH 


8o,8-7i,6 


220.078 


97 


1.25 (0.06) 


35.3 (0.2) 


8.0 (0.4) 


10.6 (0.7) 


"CO 


2-1 


220.399 


16 










CH3CN 


128-118'' 


220.476 


526 










CH3CN 


I27-II7 


220.539 


419 


0.55 (0.05) 


35.5 (0.5) 


9.7 (1.1) 


5.7 (0.9) 


HNCO 


101,9-91,8^ 


220.585 


102 


0.50 (0.04) 


33.4 (0.5) 


10.9 (1.2) 


5.8 (0.8) 


CH3CN 


126-116* 


220.594 


326 


0.84 (0.03) 


33.8 (0.3) 


15.2 (0.8) 


13.6 (0.8) 


CH3CN 


I25-II5 


220.641 


247 


0.95 (0.05) 


36.3 (0.4) 


14.1 (0.9) 


14.2 (1.3) 


C2H5CN 


252,24-242,23 


220.661 


143 


0.67 (0.03) 


35.8 (0.2) 


9.6 (0.6) 


6.8 (0.5) 


CH3CN 


I24-II4 


220.679 


183 


1.19 (0.06) 


35.2 (0.2) 


10.4 (0.6) 


13.2 (1.0) 


CH3CN 


I23-II3 


220.709 


133 


1.95 (0.08) 


35.0 (0.2) 


9.4 (0.4) 


19.5 (1.2) 


CH3CN 


I22-II2 


220.730 


97 


1.86 (0.03) 


34.9 (0.1) 


9.9 (0.2) 


19.6 (0.5) 


CH3CN 


12i-lliS 


220.743 


76 










CH3CN 


I20-II0'' 


220.747 


69 










C2H5CN 


262,25-252,24 


229.265 


154 


0.56(0.05) 


35.7 (0.4) 


8.7 (0.9) 


5.2 (0.7) 


CH3OH 


154,ii-163,l3E 


229.589 


374 


0.58 (0.05) 


35.1 (0.4) 


10.4 (1.0) 


6.5 (0.8) 


CH3OH 


8-i,8-7o,7E 


229.759 


89 


1.67 (0.07) 


35.5 (0.1) 


7.2 (0.3) 


12.8 (0.8) 


CH3OH 


195,14-204,17A- 


229.939 


579 


0.43 (0.09) 


32.1 (0.4) 


4.2 (1.0) 


1.9 (0.6) 


CH3OH 


3_2,2-4-l,4E 


230.027 


40 


0.83(0.06) 


35.8 (0.2) 


7.2 (0.6) 


6.3 (0.7) 


12CO 


2-1 


230.538 


17 










ocs 


19-18 


231.061 


111 


1.80 (0.08) 


35.0 (0.2) 


7.9 (0.4) 


15.2 (1.0) 


MM3 


CISq 


2-1 


219.560 


16 


1.54 (0.09) 


34.4 (0.1) 


3.5 (0.2) 


5.7 (0.5) 



^Formal errors from the single Gaussian fits are given in (). 
'^Blended with HNCO(102,8-92,7) 
^Blended with HNCO(102,9-92,8) 
''Profile not Gaussian-skewed to red. 

''For reference, the frequency of the k=8 CH3CN(12-11) transition is listed here and marked in Figure (2] though it is 
detected at <3cr. 

'Partially blended with nearby line. 



SBlended with CH3CN(12o-llo) 
^Blended with CH3CN(12i-lli) 
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Table 4. 1.3 mm Spectral Lines: G19.01-0.03MM1 



Fitted Line Parameters 





Tt a n c! 1 1 1 r»n 


T<Vfi(~n 1 fin r*'\7' 


F, /k 


Intensity'^ 


^ center^ 


V V IVJlill 


J uv 










Jy beam~i 


km 




km s ^ 


Jy beam km s 


C180 


2-1 


219.560 


16 


1.30 


(0.08) 


60.3 


(0.1) 


3.8 (0.3) 


5.3 (0.5) 


HNCO 


10o,i()-9o,9 


219.798 


58 


0..34 


(0.04) 


59.2 


(0.6) 


10.1 (1.5) 


3.6 (0.7) 


Hi-^CO 


3i,2-2i,i 


219.909 


33 


0..36 


(0.04) 


59.3 


(0.3) 


5.3 (0.7) 


2.0 (0.3) 


SO 


65-54'' 


219.949 


35 


0.51 


(0.06) 


59.0 


(0.2) 


3.8 (0.5) 


2.1 (0.4) 


CH3OH 


8o,8-7i,6 


220.078 


97 


0.59 


(0.06) 


59.9 


(0.3) 


5.0 (0.6) 


3.1 (0.5) 


"CO 


2-1 


220.399 


16 














CH3CN 


I24-II4 


220.679 


183 


0.34 


(0.03) 


58.5 


(0.4) 


10.0 (1.1) 


3.6 (0.5) 


CH3CN 


I23-II3'' 


220.709 


133 


0.61 


(0.05) 


60.6 


(0.3) 


6.3 (0.6) 


4.1 (0.5) 


CH3CN 


I22-II2 


220.730 


97 


0.58 


(0.05) 


60.1 


(0.2) 


3.7 (0.4) 


2.3 (0.3) 


CH3CN 


12i-lli'= 


220.743 


76 


0.58 


(0.06) 


59.5 


(0.2) 


3.9 (0.6) 


2.4 (0.4) 


CH3CN 


12o-llo'= 


220.747 


69 


0.61 


(0.05) 


59.6 


(0.2) 


4.3 (0.6) 


2.8 (0.5) 


CH30H 


8-i,8-7o,7E 


229.759 


89 


0.8 


(0.1) 


60.4 


(0.2) 


3.6 (0.5) 


3.2 (0.6) 


CH30H 


3_2,2-4-l,4Eb 


230.027 


40 


0.43 


(0.06) 


60.0 


(0.4) 


5.6 (1.0) 


2.56 (0.6) 


12CO 


2-1 


230.538 


17 














OCS 


19-18 


231.061 


111 


0.60 


(0.06) 


60.0 


(0.3) 


5.0 (0.6) 


3.2 (0.5) 



'^Formal errors from the single Gaussian fits are given in (). 

''Profile not Gaussian. CH3OH lines have blue wings. CH3CN k=3 line is skewed to blue of fit. 

■^CHsCN (I20-II0) and CH3CN (12i-lli) are partially blended. Parameters for CH3CN k=0 and k=l lines 
derived from a two-component Gaussian fit. 



Table 5. CH3CN Fit Parameters 



Source Distance 


T 




Size 




N(CH3CN) 


Abundance 


N(H2) 


n(H2) 






kpc 


K 




pc 


AU 


cm~^ 






cm~^ 


cm~3 


M0 


Single Component Fit 
























G11.92-0.61-MM1 


3.8 


140 


1.0 


0.02 


3800 


3.3x10^8 


10" 


7 


3.3x10^3 


1.2x10^" 


1.9 
















10- 


8 


3.3x1024 


1.2x10^ 


19 
















10" 


9 


3.3x10^5 


1.2x10^ 


192 


G19.01-0.03-MM1 


4.2 


114 


0.6 


0.01 


2500 


1.7x10^^ 


10" 


7 


1.7x10^3 


8.8x10*^ 


0.4 
















10" 


8 


1.7x1024 


8.8x10'^ 


4.2 
















10" 


9 


1.7x10^6 


8.8x10* 


42 


Two Component Fit 
























G11.92-0.61-MM1 
























Cool 


3.8 


77 


3.0 


0.06 


11400 


l.lxlOi'S 


10- 


7 


1.1X10^2 


1.3 xlO^ 


0.6 
















10" 


8 


1.1X1023 


1.3x10^ 


5.6 
















10" 


9 


l.lXl024 


1.3x10'^ 


56 


Warm 


3.8 


166 


0.6 


0.01 


2300 


2.3x10" 


10- 


7 


2.3x1024 


1.4x10* 


4.9 
















10" 


8 


2.3X102'"'' 


l.dxlO'' 


49 
















10" 


9 


2.3xl02S 


1.4xl0i" 


490 
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Table 6. Derived Properties of Millimeter Continuum Sources 



Source 




'^dust 


Mgas 










K 




M0 


1023 cm-2 


10'' cm-3 




SMA 1.3 mm 


Gil. 92 


-0.61-MMl 


70-190 


0.026-0.010 


23.5-8.2 


6.6-2.3 


1.0-0.3 


Gil. 92 


-0.61-MM2 


20-40 


0.056-0.028 


61.6-26.3 


17.4-7.4 


2.6-1.1 


G11.92- 


-0.61-MM3 


30-80 


0.009-0.003 


8.8-2.9 


2.5-0.8 


0.4-0.1 


Total 








93.9-37.4 






5a 




20/50 




5.2/1.7 






G19.01- 


-0.03-MMl 


100-130 


0.016-0.012 


15.9-12.0 


3.9-3.0 


0.5-0.4 


5a 




20/50 




6.3/2.1 












CAKAIA 1.1 


miii 






Gil. 92 


-0.61-MMl 


70-190 


0.066-0.024 


12.9-4.5 


16.8-5.8 


5.3-1.8 


G11.92- 


-0.61-MM2 


20-40 


0.175-0.084 


41.6-17.3 


53.9-22.4 


16.9-7.1 


G11.92- 


-0.61-MM3 


30-80 


0.033-0.012 


7.0-2.3 


9.1-3.0 


2.9-1.0 


Total 








61.5-24.1 






5a 




20/50 




6.9/2.3 












CARMA 3.4 


mm 






G19.01- 


-0.03- MMl 


100-130 


0.002-0.001 


40.2-30.7 


1.9-1.5 


0.1-0.09 


5a 




20/50 




20.6/7.7 







'■Beam-averaged quantities 
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Table 7. Outflow Parameters 



Vmin^ ^max'^ Min( (l)- Vigfl [ ) Mout 0^° Pout Eout Length Tdyn Mout Pout 

km km s"'^ km ° km ergs pc yrs yr~^ km yr~^ 

xlO^B xlO-4 



Gll.92-0.61: Main Outflow 



CO(2-l) 



Red 48.1 


71.2 


13 


0.2 


- 


3.1 


0.6 


0.2 


5000 


0.4 


0.001 










10 


17.7 


19.0 


0.2 


900 


2.0 


0.020 










30 


6.1 


2.3 


0.2 


2900 


0.6 


0.002 










60 


3.5 


0.8 


0.3 


8700 


0.2 


0.0004 


Blue -24.4 


21.8 


13 


0.6 




12.6 


3.2 


0.4 


6000 


1.0 


0.002 










10 


72.6 


107 


0.4 


1100 


5.6 


0.069 










30 


25.2 


12.9 


0.4 


3500 


1.7 


0.007 










60 


14.6 


0.4 


0.8 


10400 


0.6 


0.001 


Total"^ 






0.8 




15.7 


3.8 


0.5 




1.3 


0.003 










10 


90.3 


126 


0.6 




7.6 


0.089 










30 


31.4 


15.2 


0.6 




2.3 


0.009 










60 


18.1 


5.1 


1.1 




0.8 


0.002 


HCO+(1-0) 






















Red 48.1 


71.1 


13 


2.6 




55.2 


12.9 


0.2 


5000<^ 


5.2 


0.011 










10 


317.7 


427 


0.2 


900 


29.5 


0.360 










30 


110.3 


51.5 


0.2 


2900 


9.0 


0.038 










60 


63.7 


17.2 


0.4 


8700 


3.0 


0.007 


Blue -14.3 


21.8 


13 


5.2 




134.2 


40.2 


0.4 


6000<i 


8.7 


0.022 










10 


773.0 


1330 


0.4 


1100 


49.4 


0.731 










30 


268.5 


161 


0.4 


3500 


15.1 


0.078 










60 


155.0 


53.6 


0.7 


10400 


5.0 


0.015 


Total<= 






7.8 




189.4 


53.1 


0.6 




13.9 


0.033 










10 


1090.7 


1760 


0.6 




78.9 


1.091 










30 


378.8 


212 


0.6 




24.1 


0.116 










60 


218.7 


70.8 


1.1 




8.0 


0.022 



Gll.92-0.61: Northern Outflow 



CO(2-l) 



Red 


48.1 64.6 13 


0.14 




2.9 


0.6 


0.4 


12000 


0.1 


0.0002 








10 


16.6 


20.7 


0.4 


2100 


0.7 


0.008 








30 


5.8 


2.5 


0.4 


6900 


0.2 


0.001 








60 


3.3 


0.8 


0.7 


20800 


0.07 


0.0002 


Blue 


8.6 21.8 13 


0.05 




0.8 


0.1 


0.1 


5000 


0.09 


0.0002 








10 


4.6 


4.8 


0.1 


900 


0.5 


0.005 








30 


1.6 


0.6 


0.2 


2900 


0.2 


0.001 








60 


0.9 


0.2 


0.3 


8700 


0.05 


0.0001 


Total 




0.2 




3.7 


0.7 


0.5 


8500== 


0.2 


0.0004 








10 


21.3 


25.5 


0.5 


1500 


1.3 


0.014 



Table 7— Continued 



Min{\v—VLSR\) Mout Pout Eout Length Tj_y„ Mout Pout 

km km km s"-' Mq ° Mq km ergs pc yrs Mq yr^^ Mq km s^^ yr~^ 



30 7.4 3.1 0.6 4900 0.4 0.002 

60 4.3 1.0 1.0 14700 0.1 0.0003 

HCO+(1-0) 



Red 


48.1 


59.6 13 


0.8 




13.4 


2.4 


12000=^ 


0.7 


0.001 










10 


77.1 


79.3 


2100 


3.7 


0.036 










30 


26.8 


9.6 


6900 


1.1 


0.004 










60 


15.5 


3.2 


20800 


0.4 


0.001 


Blue 


15.3 


21.8 13 


0.3 




4.0 


0.6 


5000<i 


0.5 


0.001 










10 


23.1 


21.2 


900 


2.9 


0.026 










30 


8.0 


2.6 


2900 


0.9 


0.003 










60 


4.6 


0.9 


8700 


0.3 


0.001 


Total 






1.0 




17.4 


3.0 


8500=^ 


1.2 


0.002 










10 


110.2 


100 


1500 


6.9 


0.067 










30 


34.8 


12.1 


4900 


2.1 


0.007 










60 


20.1 


4.0 


14700 


0.7 


0.001 



G19.01-0.03 



CO(2-l) 



Red 


75.6 


88.8 


16 


0.1 




2.3 


0.5 


0.4 


16400 


0.07 


0.0001 












10 


13.2 


15.3 


0.4 


2900 


0.41 


0.005 












30 


4.6 


1.9 


0.5 


9500 


0.13 


0.0005 












60 


2.6 


0.6 


0.9 


28400 


0.004 


0.00009 


Bluef 


-46.4 


39.3 


21 


0.3 




8.4 


3.3 


0.3 


3000 


1.0 


0.003 












10 


48.6 


110 


0.3 


500 


6.0 


0.09 












30 


16.9 


13.2 


0.4 


1700 


1.8 


0.01 












60 


9.8 


4.4 


0.6 


5200 


0.6 


0.002 


NE Blue clump 


9.6 


39.3 


21 


0.1 




2.2 


0.7 










Totals 








0.5 




12.9 


4.5 


0.7 


5000 


1.1 


0.003 












10 


74.3 


149 


0.7 


900 


6.4 


0.08 












30 


25.8 


17.9 


0.8 


2900 


2.0 


0.009 












60 


14.9 


6.0 


1.4 


8700 


0.7 


0.002 


HCO+(1-0) 
























Red 


66.9 


83.4 


6.9 


3.6 




40.4 


5.4 


0.4 


16400^^ 


2.2 


0.002 












10 


232.4 


177 


0.4 


2900 


12.5 


0.080 












30 


80.7 


21.4 


0.4 


9500 


3.8 


0.009 












60 


46.6 


7.1 


0.7 


28400 


1.3 


0.002 


Blue 


7.9 


53.8 


6.2 


9.5 




155.7 


42.3 


0.4 


3000<^ 


31.5 


0.052 












10 


896.4 


1400 


0.4 


500 


179 


1.695 












30 


311.3 


169 


0.5 


1700 


54.5 


0.180 












60 


179.7 


56.4 


0.8 


5200 


18.2 


0.035 


Total 








13.1 




196.0 


47.7 


0.8 


5000 


26.1 


0.039 












10 


1128.7 


1580 


0.8 


900 


148 


1.280 
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Table 7 — Continued 





Vmin'' Vmax'' M'mQv-VLS r\) 






Pout 


Eout 


Length 


^dyn 


Mout 


Pout 




km s km s km s 


Mq 




Mq km 


ergs 


pc 


yrs 


M© yr-i 


Mq km yr"-^ 












X lO''-' 
















30 


392.0 


191 


0.9 


2900 


45.2 


0.136 








uu 
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uo.o 


1 ^ 


o / uu 




n n9R 
u.uzu 


rtea 


Oo.D 00.4 O.O 






Q1 Q 

oi.y 


4.0 






1.0 


n nn9 
U.UUz 








1 n 


1 8*^ 7 
ioo- / 


iOo 






o.o 










'in 


Do.o 


1 Q 1 






9 ^ 
Z.O 


u.uu / 








OU 


OD.O 


R 4 






U.O 


U.UUl 






D.O 












91 ^ 
Zl.O 


n OAR 








1 n 


7Qn 1 








1 99 
J.ZZ 


1 4Qzl 








30 


274.4 


165 






37.3 


0.158 








60 


158.4 


54.9 






12.4 


0.030 


Total 




8.9 




169.1 


45.9 






17.7 


0.034 








10 


973.8 


1520 






100 


1.104 








30 


338.2 


184 






30.7 


0.117 








60 


195.2 


61.3 






10.2 


0.023 


Red 


70.2 83.4 10.2 


1.8 




26.4 


4.3 






1.1 


0.002 








10 


151.8 


142 






6.0 


0.052 








30 


52.7 


17.2 






1.8 


0.006 








60 


30.4 


5.7 






0.6 


0.001 


Blue 


7.9 50.5 9.5 


5.0 




125.6 


40.3 






16.6 


0.042 








10 


723.4 


1340 






94.2 


1.367 








30 


251.2 


161 






28.8 


0.145 








60 


145.0 


53.7 






9.6 


0.028 


Total 




6.7 




152.0 


44.6 






13.5 


0.030 








10 


875.2 


1480 






76.3 


0.993 








30 


303.9 


178 






23.3 


0.105 








60 


175.5 


59.4 






7.8 


0.020 



] 



'^Central velocity of first/last channel used for calculation of outflow properties. 
''Inclination angle from the plane of the sky. 
"^Sum of red and blue lobes, 
'iprom I^CO 

''Calculated from the end-to-end outflow length divided by 2 and 'Vmax,outflow=^^ km s~^. 

^Excluding channels v=26.1 and 22.8 km s~'. These channels are dominated by image ajtifax;ts from resolved-out 
large scale emission. 

8 Including NE blue clump 



Right Ascension (J2000} 



Right Ascension (J2CKK)) 



Fig. 1. — Three color Spitzer images (3.6, 4.5, 8.0 /xm: blue, green, red) of the EGOs (a) 
Gil. 92— 0.61 and (c) G19.01— 0.03 overlaid with contours of 1.3 mm continuum emission 
from the SMA (black contours, levels 5,10,30 x a = 3.5 mJy beam~^). The SMA beam 
is shown at lower left in each panel. In all panels, red diamonds mark positions of 6.7 
GHz CH.sOH mas e rs an d magenta crosses mark positions of 44 GHz CH3OH masers from 



Cyganowski et al.l (120091 ). The 8 /xm nebula to the NE of Gil. 92— 0.61 (extreme upper left in 



panel (a)) is associated with IRAS 18110-1854. The black rectangle in (a) indicates the field 
of view of (b). (b) The inset shows the 24 yum emission towards Gil. 92— 0.61 in greyscale, 
overlaid with contours of 1.4 mm continuum emission (blue) from CARMA (levels 4,5,10,20 
X cr = 4.3 mJy beam~^) and 8 fim emission (red) from GLIMPSE (levels 250,500 MJy 
sr-i). G11.92-0.61-MM1, MM2, and MM3 are labeled "1", "2", and "3", respectively. The 
resolution of the 24 fim image is 6" (sc alebar) and the CARMA bea m is shown at lower left. 
The position of the H2O maser from IHofner fc Churchwelll (119961 ) is marked with a filled 
green diamond. The 24 /im emission towards Gil. 92— 0.61 is saturated in the MIPS image 
(white pixels). In (c), the MIPS 24 /im emission towards G19.01— 0.03 is shown as yellow 
contours (levels 900, 1500 MJy sr'^). 
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Fig. 2. — SMA LSB and USB spectra towards the 1.3 mm continuum peak of 
G11.92— 0.61MM1. The spectra have been Hanning smoothed. Lines detected at >3cr are 
labeled and listed in Table |3l 
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Right Ascension (^2000) 

Fig. 3. — Gil. 92— 0.61: The greyscale shows the SMA 1.3 mm continuum emission. Contours 
are drawn from integrated intensity (moment zero) maps of the indicated molecular line. The 
species and upper state energy in K of the transition are listed at upper right in each panel. 
Contour levels are: C^^O: 1.92, 3.84, 5.76, 7.68 Jy bcam-^ km s'^ H^^CO: 1.65, 3.3, 4.95 
Jy beam-i km s^^ SO: 3.3, 4.95, 13.2, 23.1 Jy bcam^^ km s^^ HNCO: 2.76, 5.52, 8.28, 
11.04 Jy beam-i km s"^; CH3CN (k=2): 2.82, 8.46, 14.1 Jy beam'^ km s^^; OCS: 2.4, 7.2, 
12 Jy beam-i km s^^; C2H5CN: 1.62, 3.24, 4.86 Jy beam"^ km s'^; CH3CN (k=7): 1.8, 
3.6 Jy beam~^ km s~^; HC3N(v7=l): 2.52, 5.04 Jy beam~^ km s~^. These contour levels 
are (3,6,9,12) xa for C^^O and HNCO; (3,6,9) xa for H^^CO and C2H5CN; (3,9,15) xa for 
CH3CN (k=2) and OCS; (3,6) xa for CH3CN (k=7) and HC3N; and (3,4.5,12,21) xa for 
SO. The SMA beam is shown at lower left. 
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Right Ascension (J2000) 



20 40 _ 60 

LSR Velocity (km s ) 



Fig. 4. — Gil. 92— 0.61: (a) The greyscale shows the SMA 1.3 mm continuum emission. 
Contours are drawn from integrated intensity (moment zero) maps of the indicated CH3OH 
hne. The upper state energy in K of the transition is given at upper right in each panel. 



Contour levels are: E, 
10.08 Jy beam 



upper 

^ km s~^; E, 



40 K: 2.2, 4.4 Jy beam"^ km s"^; E, 



upper " 



upper' 



--89K: 2.52, 5.04, 7.56, 
=97 K: 2,4,8 Jy beam'^ km s'^ E„pper=374 K: 2.8, 5.6 Jy 

upper 40 K and ^upper 374 K 



beam~^ km s~^. These contour levels are (4,8) xa for the E, 
transitions, (4,8,12,16) xa for the E„pper=89 K transition, and (4,8,16) xa for the Eupper=97 
K transition. Black cro sses mark the positions of 44 GHz Class I CH3OH masers from 
Cyganowski et al.l ( 120091 ). The SMA beam is shown at lower left, (b) and (c): Spectra at the 
peak positions of the indicated CH3OH (8_i,8-7o,7) emission. Solid line: 8-1^-70,7 transition 
(229.759 GHz, E„pper=89 K); dashed line: 3_2,2-4_i,4 (230.027 GHz, E„ppe^=40 K); dotted 
line: 8o,8-7i,6 (220.078 GHz, E„ppe^=97 K). 
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Fig. 5. — Gil. 92— 0.61: First moment maps of the indicated species (colorscale) overlaid with 
contours of the SMA 1.3 mm continuum emission (levels 10, 30 x a = 3.5 mJy beam~^). 
The inner contour is approximately the size of the SMA beam. E^pper in K is listed under the 
molecule nam e in the upper right o f each panel. The open triangle marks the water maser 
position from lHofner fc Churchwelll f|l996r). The filled squa re marks the position of the 6.7 
GHz Class II CH3OH maser from ICyganowski et al.l (120091 ). 
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Fig. 6. — Gil. 92-0. 61: Three color Spitzer images (3.6, 4.5, 8.0 /im: blue, green, red) overlaid 
with contours of 1.4 mm continuum emission (black) and high velocity (a) ^^C0(2-l), (b) 
HCO"'"(1-0), and (c) SiO(2-l) emission. In each panel, positions of 6.7 GHz CH3OH masers 
are marked with black diamonds, and po sitions of 44 GHz CH3OH masers are marked with 
magneta crosses (jCyganowski et ah 2009 ). Continuum contour levels are 4,5,10,20 x a = 4.3 
mJy beam~^. The vlsr is ~35 km s~^ ( §3.1.3|l4.2p . (a) ^^C0(2-l) emission integrated over 
v=-24.4 to 21.8 km s"^ (blue) and v=48.1 to 71.2 km s~^ (red), e.g. ~t;L5ij±13 km s"^ 
Levels are 5,9,17,29,45 Jy beam~^ km s~^ for both red and blue contours. The SMA and 
CARMA beams are shown at lower left, (b) HCO"*"(1-0) emission integrated over v=-14.3 
to 21.8 km s~^ (blue) and 48.1 to 71.1 km s~^ (red), e.g. ~fLS'R±13 km s~^. Contour levels: 
Blue: 1.3, 1.8, 2.8, 4.3, 5.8 Jy beam'^ km s'^; Red: 1.3, 1.8, 2.8, 3.8 Jy beam'^ km s-\ 
The CARMA beams are shown at lower left, (c) SiO(2-l) emission integrated over v=8.0 
to 24.9 km s~^ (blue) and v=45.1 to 53.5 km s~^ (red), e.g. ~fL5'/j±10 km s~^. Contour 
levels: Blue: 0.8, 1.0, 1.4, 2.0 Jy beam'^ km s'^ Red: 0.6, 0.8, 1.0 Jy beam-^ km s'K The 
CARMA beams are shown at lower left. 
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Fig. 7.— GU. 92-0.61: Channel maps of ^^00(2-1). The greyscale is the 1.3 mm SMA 
continuum image and the contours are the ^^CO emission (levels 0.24, 0.48, 0.96, 1.44, 2.16, 
3.6 Jy beam~^ (solid), —0.24, —0.48 Jy beam~^ (dotted)). Each panel is the average of two 
channels of the smoothed ^^CO cube; the center velocity of each panel is given at upper 
right. The contour levels are (—8,-4,4,8,16,24,36,60) x the rms in a line-free panel. The 
field of view shown is the same as in Figures [8] and [91 The dashed circle shows the FWHP 
of the SMA primary beam. The SMA synthesized beam is shown as a filled ellipse at lower 
left. 
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Fig. 8.— Gll.92-0.61: Channel maps of HCO+(1-0). The greyscale is the 1.3 mm SMA 
continuum image. The contours are the HCO+ emission (levels 0.06, 0.096, 0.144, 0.192, 
0.288, 0.432, 0.576, 0.864 Jy beam-^ (solid), -0.06, -0.096 (dotted)). The contour levels 
are (-8,-5,5,8,12,16,24,36,48,72) x the rms in a line-free panel. The center velocity of each 
panel is given at upper right. The field of view shown is the FWHP primary beam of the 
10 m CARMA antennas (80" square). The CARMA synthesized beam is shown as a filled 
ellipse at lower left. 
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Fig. 9.— Gll.92-0.61: Channel maps of SiO(2-l). The greyscale is the 1.3 mm SMA contin- 
uum image. The contours are the SiO emission (levels 0.055, 0.088, 0.132, 0.176, 0.264, 0.396 
Jy beam^-*^ (solid), —0.055 Jy beam" -"^ (dotted)). The contour levels are (—5,5,8,12,16,24,36) 
X the rms in a line-free panel. The center velocity of each panel is given at upper right. The 
field of view shown is the FWHP primary beam of the 10 m CARMA antennas (80" square). 
The CARMA synthesized beam is shown as a filled ellipse at lower left. 
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Fig. 10. — SMA LSB and USB spectra towards the 1.3 mm continuum peak of 
G19.01— 0.03MM1. The spectra have been Hanning smoothed. Lines detected at >3cr are 
labeled and listed in Table HI 
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Fig. 11. — G19.01— 0.03: The greyscale shows the SMA 1.3 mm continuum emission. Con- 
tours are drawn from integrated intensity (moment zero) maps of the indicated molecular 
line. The species and upper state energy in K of the transition are listed at upper right in 
each panel. Contour levels are: C^^O: 1.65, 3.3, 4.95 Jy beam'^ km s^^; H^^CO: 1.2, 2.0 
Jy beam"^ km s"^; SO: 0.9, 1.5 Jy beam"^ km s"^; HNCO: 1.32, 2.20 Jy beam"^ km s"^; 
CH3CN (k=2): 1.11, 1.85 Jy beam-^ km 8"^ OCS: 1.2, 2.4 Jy beam'^ km 8"^ The8e 
contour level8 are (3,6,9) x a for C^^O, (3,5) x a for H^^CO, SO, HNCO, and CH3CN, and 
(3,6) X (7 for OCS. The SMA beam i8 8hown at lower left. 




Fig. 12. — G19.01— 0.03: The greyscale shows the SMA 1.3 mm continuum emission. Con- 
tours are drawn from integrated intensity (moment zero) maps of the indicated CH3OH hne. 
The upper state energy in K of the transition is given at upper right in each panel. Con- 
tour levels are: E„pper-=40 K: 1.11, 1.85 Jy beam~^ km s~^; E„pper=89 K: 1.02, 1.7, 2.38 Jy 
beam~^ km s^^; Eupper=^7 K: 1.32, 2.20, 3.08 Jy beam~^ km s~^. These contour levels are 
(3,5) X a for the Eupper=40 K transition and (3,5,7) x a for the E„pper=89 K and E„ppe.r=97 
K transitions. Bl a ck cro sses mark the positions of 44 GHz Class I CH3OH masers from 
Cyganowski et al.l ( 120091 ). The SMA beam is shown at lower left. Spectra are shown at 
the peak position of the indicated CH3OH (8_i,8-7o,7) emission spot. Solid line: 8_i,8-7o,7 
transition (229.759 GHz, E„pper=89 K); dashed line: 3_2,2-4-i,4 (230.027 GHz, E„ppe^=40 K); 
dotted line: 8o,8-7i,6 (220.078 GHz, E„ppe^=97 K). 
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Fig. 13. — G19.01-0.03: Three color Spitzer images (3.6, 4.5, 8.0 /im: blue, green, red) over- 
laid with contours of high velocity (a) ^^00(2-1), (b) HCO+(1-0), and (c) SiO(2-l) emission 
In each panel, positions of 6.7 GHz CH3OH masers are marked with bl ack diamonds, and 



posit ions of 44 GHz CH3OH masers are marked with magneta crosses (jCyganowski et al. 



20091). The vlsr is ~60 km s~^ (gSXS]). (a) 12C0(2-1) emission integrated over v=-46.4 



to 19.5 km s~^ (blue) and v=75.6 to 88.8 km s~^ (red). Contour levels: Blue: 7.2, 9.6, 12, 
15.6, 19.2, 22.8 Jy beam"^ km s~^; Red: 4.8, 7.2, 9.6 Jy beam~^ km s~^. Black contours: 
SMA 1.3 mm continuum emission (levels 5,10,30 x cr = 3.5 mJy beam~^). The SMA beam 
is shown at lower left, (b) HCO"''(1-0) emission integrated over v=7.8 to 52.2 km s~^ (blue) 
and 68.6 to 83.4 km s~^ (red), e.g. ~fi5R±8 km s~^. Contour levels: Blue: 1.0, 1.2, 1.6, 2.2, 
3.0, 4.0 Jy beam"^ km s"^; Red: 0.6, 0.8, 1.0, 1.2 Jy beam'^ km s'^ The CARMA beam is 
shown at lower left, (c) SiO(2-l) emission integrated over v=26.0 to 52.9 km s~^ (blue) and 
v=66.4 to 83.3 km s~^ (red), e.g. ~fL5_R±7 km s~^. Contour levels: Blue: 0.5, 0.8, 1.2, 1.8, 
2.6 Jy beam^^ km s^^; Red: 0.5 Jy beam~^ km s~^. The CARMA beam is shown at lower 
left. 
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Fig. 14.— G19.01-0.03: Channel maps of ^2CO(2-l). The greyscale is the 1.3 mm SMA 
continuum image and the contours are the ^^CO emission (levels 0.2, 0.4, 0.6, 0.8, 1.2, 1.8, 
3 Jy beam~^(sohd), —0.2, —0.4 Jy beam~^ (dotted)). Each panel is the average of three 
channels of the smoothed ^^CO cube; the center velocity of each panel is given at upper 
right. The contour levels are (-8,-4,4,8,16,24,36,60) x the rms in a line-free panel. The 
field of view shown is the same as in Figures [15] and [161 The dashed circle shows the FWHP 
of the SMA primary beam. The SMA synthesized beam is shown as a filled ellipse at lower 
left. 
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Fig. 15.— G19.01-0.03: Channel maps of HCO+(1-0). The greyscale is the 1.3 mm SMA 
continuum image. The contours are the HCO^ emission (levels 0.03, 0.048, 0.072, 0.096, 
0.144, 0.216, 0.288, 0.432, 0.648, 0.864 Jy beam-^ (solid), -0.03, -0.048 Jy beam-i(dotted)). 
The contour levels are (-8,-5,5,8,12,16,24,36,48,72,108,144) x the rms in a line-free panel. 
The center velocity of each panel is given at upper right. The velocity range shown matches 
that of Figure [TH Note that the velocity coverage of the CARMA HCO^ observations does 
not extend blueward of ~8 km s~^. The field of view shown is the FWHP primary beam 
of the 10 m CARMA antennas (80" square). The CARMA synthesized beam is shown as a 
filled ellipse at lower left. 
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Fig. 16. — G19.01— 0.03: Channel maps of SiO(2-l). The greyscale is the 1.3 mm SMA 
continuum image. The contours are the SiO emission (levels 0.03, 0.048, 0.072, 0.096, 
0.12, 0.144,0.216,0.288 Jy beam"! (solid), -0.03 Jy beam^^ (dotted)). The contour levels 
are (—5,5,8,12,16,20,24,36,48) x the rms in a line-free panel. The center velocity of each 
panel is given at upper right. The velocity range shown matches that of Figure [TH Note 
that the velocity coverage of the CARMA SiO observations does not extend blueward of ~8 
km s~^. The field of view shown is the FWHP primary beam of the 10 m CARMA antennas 
(80" square). The CARMA synthesized beam is shown as a filled ellipse at lower left. 
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Fig. 17. — Gil. 92— 0.61: (a) Data (black circ les) overlaid with the b est fit model (solid black 
line) and other good fits (grey lines) from the lRobitaille et al.l ( 120071 ) model grid. The twenty 
best fits are shown, (b) Histogram of the bolometric luminosity of the 20 best-fit models 
(hatched histogram) and all models in the grid (grey histogram). 
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Fig. 18. — G19.01— 0.03: (a) Data (black circles) overlaid with the bes t fit model (solid 
black line) and other good fits (grey lines) from the iRobitaille et al.l (j2007| ) model grid. The 
twenty best fits are shown. The dashed line shows the stellar photosphere, in the absence 
of circumstellar dust, of the central source for the best-fit model. Histograms of (b) the 
bolometric luminosity, (c) the envelope mass, and (d) the stellar age for the 20 best-fit 
models (hatched histograms) and all models in the grid (grey histograms). 
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Fig. 19. — G19.01— 0.03: Envelope accretion rate v. stellar mass for the 20 best- fit models 
(black dots) and all models in the grid (grey dots) . 
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Fig. 20.— SMA CH3CN(12-11) spectrum towards the G11.92-0.61-MM1 continuum peak 
(black histogram) overlaid with the best-fit one component (red line) and two component 
(blue line) models. The physical parameters of the best-fit models are summarized in 
Table [51 For the single component model, ^mfwhm=Q km s~^; for the two component 
model, Avfvkhm(coo1)=7 km s~^and /S.Y pwHAiiyf 8iim)=Q km s^^. For all components, 
^LSR=35.5 km s~^. The CH3CN model accounts for emission from CHg^CN, assuming 
CH3^C:CH3^C=60:1. The strongest observed CHg^CN lines, including those blended with 
the CH3CN k=5 and k=6 components, are indicated. The C2H5CN line at 220.661 GHz is 
not fit, nor is the HNCO line at 220.585 GHz (partially blended CH3CN k=6). 
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Fig. 21.— SMA CH3CN(12-11) spectrum towards the G19.01-0.03-MM1 continuum peak 
(black histogram) overlaid with the best-fit single component model ( red line). The physical 
parameters of the best-fit model are summarized in Table O 
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Fig. 22. — Rotation diagrams for the CH3OH transitions observed at the G11.92— 0.61-MMl 
continuum peak. The upper panel shows the fit to the raw data. In the lower panel, the 
column densities have been corrected for optical depth effects as described in §4.1.21 The 
fitted temperatures are indicated in each panel. 



